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ABSTRACT 

We present the rest-frame optical spectral properties of 155 luminous quasars at 3.3 < z < 6.4 taken 
with the AKARI space telescope, including the first detection of Ho emission line as far out as 
z ^ 6. We extend the scaling relation between the rest-frame optical continuum and line luminosity of 
active galactic nuclei (AGNs) to the high luminosity, high redshift regime that has rarely been probed 
before. Remarkably, we find that a single log-linear relation can be applied to the 5100 A and Ha 
AGN luminosities over a wide range of luminosity (10'^^ < L 5100 < 10'^^ ergs s“^) or redshift (0 < z < 6), 
suggesting that the physical mechanism governing this relation is unchanged from z = 0 to 6, over five 
decades in luminosity. Similar scaling relations are found between the optical and the UV continuum 
luminosities or line widths. Applying the scaling relations to the H/3 black hole mass (Mbh) estimator 
of local AGNs, we derive the Mbh estimators based on Ha, Mg 11, and G IV lines, finding that the 
UV-line based masses are overall consistent with the Balmer-line based, but with a large intrinsic 
scatter of 0.40 dex for the C IV estimates. Our 43 Mbh estimates from Ha confirm the existence of 
BHs as massive as ~ IO^^Mq out to z ~ 5, and provide a secure footing for previous Mg 11-line based 
studies that a rapid Mbh growth has occurred in the early universe. 

Subject headings: galaxies: active — galaxies: evolution — quasars: emission lines — quasars: super- 
massive black holes 


1. INTRODUCTION 


Quasars, galaxies that are in an active phase due to vig¬ 
orous accretion of matter toward the central supermas- 
sive black hole ( BH), has been vast l y disc o vered by many 


survevs fe.g.. ISchmidt & Greeiil 

1983 

: iHewett et al. 

1995 

IBovle et al. 

120001: lYork et al.l 

2000 

[Richards et al. 

2002 

Ilm et al.l 

20071 iLee et al.l 

2008 

; iWillott et al. 


2010bl : iWu et alJ 1201011 . Through the discovery of 
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quasars 

at high redshift dPan et al.1 120001: ICool et al. 

2006 

; 

Gotol 120061: IStern et a 

.1 120071: IWillott et al. 

2007 

Mortlock et al.l 120111: 

Venemans et al.l 120131: 

Banados et al.ll2014f). we are witnessing the early stages 


of supermassive BH growth in the distant universe. The 
number density of optically luminous quasars at high 
redshift, quickly increases w ith cosmic time towards 
its maximum at z = 2'^3 (iDunlop fc Peacock! 119901: 
Warren et al.lll994 iSchmidt et al.lll995HKennefick et al.l 

19951 : [Richards et al.l l2006bl : iMcGreer et al.l 1201311 . 
Accompanied by high accretion rates amon g lumi¬ 
nous quasars at z>A fe.g.. FWillott et al.l l2010at 
iDe Rosa et al.ll20lll iDe Rosa et al.1120141 1. this suggests 
a rapid BH growth of the luminous population of ac¬ 
tive galactic nuclei (AGNs) in the early universe. Also, 
an unusual population of AGNs known as dust-poor 
quasars - quasars with little infrared emission from hot 
and war m dust - is found to be more common at higher 
redsh ift (j.liang- et al.ll20inl : l.lun fc Imll2013l : [LeiDski et al.l 
l2014f l. Luminous dust-poor quasars tend to have lower 
BH masses (Mbh) or higher Eddington ratios compared 
to typ ical luminous quasars (|Jiang et al.ll201Cll : lJun fc Iml 
l2013l) . indicating the build-up of the AGN dusty sub¬ 
structure during its early mass accretion. 


One of the key findings in the study of high red¬ 
shift quasars is that there exist extremely massive BHs, 
with the mass reac hing Mrh ^ 10^^Mp) at z = 2-5, and 


^ lO^Mrn at z =6-7 (IJiang et al.ll2007:INetzer et al.ll2007l: 

Kurk et al.l 120071: IVester^aard et al. 2008 

: Shen et al.1 

20081: IMortlock et al.l 120111: IDe Rosa et al. 

1201411. IJn- 


der the concordance cosmology, the time gap between the 
reionization epoch of the universe from the recent Planck 
study of the cosmic microwave b ackground, z ^ 11.5 
(|Planck Collaboration et ^l2014ll . and z = 6 is 0.5 Gyr. 
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Considering the case where a Population II stellar seed 
Mbh starts to grow at z ~ 11 .5, the given time is too 
short for the seed to become an extremely massive BH 
at z = 6. Under the Eddington-limited accretion where 
the mass accretes at a maximal rate with the radiative 
pressure gradient and gravity in balance, we expect a BH 
to grow as 


M{t) = Mo (I) 


where Mq is the mass at toj e is the radiative efficiency, 
and tEdd is the Eddington limited timescale of 0.45 Gyr. 
Assuming a typical value of e = 0.1, a BH can grow by 
2x10^ times over the time span of 0.5 Gyr, without 
considering feedback mechanisms that could slow down 
the BH growth. This maximal growth factor is far too 
small for a stellar mass BH with a typical seed mass of 
IOM 0 , to grow into the extremely massive AGNs that 
have been observed recently. Consequently, BH seeds 
that may have started accreting prior to the reionization 
epoch j which are more massive (e.g.. iBromm &: Lo^ 
200.‘U iBegelman et al.ll2005 iLodato fc Nataraianll2006F 


Bellovarv et al.l 201 ID or go through super-Eddington 


accret ion fe.g.. IVolonteri fc Reesll2005t iWvithe fc LoebI 
[Ml IMadau et al.ll20I4ll . are suggested to explain the 
Mbh of quasa rs at high redshi ft (also see a review on 
this subject bv lNataraiaDll20l3l . 

Obviously, accurate determination of Mbh is an im¬ 
portant requirement for understanding the BH growth 
at high redshift. This is especially true for BHs at the 
most massive end. A large uncertainty in Mbh can scat¬ 
ter the abundant lower mass BHs into the high mass end 
of Mbh distribution, while the effect in the opposite di¬ 
rection is much less significant since higher mass BHs are 
relatively rare. As a result, the number of extremely mas¬ 
sive BHs can be easily overestimated. In principle, the 
correction to this effect is possible, but it requires a good 
knowledge on the error of Mbh measurements which is, 
however, rather difficult to obtain. This poses a poten¬ 
tial challenge to the understanding of the BH growth at 
high redshift as described below. 

In most of high redshift quasar studies, BH masses 
are estimated using the rest-frame ultraviolet (UV)-part 
of spectra which is redshifted into the rest-frame op¬ 
tical (rest-optical). The velocity widths of broad UV 
lines such as G IV and Mg II are used as measures 
of the gas motion of the broad line region (BLR), and 
the continuum or the line luminosity at the rest-frame 
UV (rest-UV) is used as a proxy for the size of BLR 
(.Rblr)- One gets Mbh by combining the two pieces 
of information through a virial mass estimator, Mbh oc 
Rr t.r X EWHMbt p (e.g.. iMcLure fc Dunlod 1^004 


Vestergaard 1 12004 iBaskin fc Laor 20051: iSulentic et SI 

2007t ISben et aLM200Sl IPa.rk et a1.l 12071 . While UV- 


based Mbh estimators are useful tools to measure the 
Mbh of AGNs, they are secondary estimators that are de¬ 
rived from rest-optical spectral properties such as the line 
luminosity or width of H/3 and Ha, and the continuum 
luminosity at 5100 A (Laioo)- Assuming the UV luminos¬ 
ity follows the optical broad line region radius-luminosity 
{Rblr-L) relation with a constant factor, and the UV 
broad line width follows the optical line width as a power 
law relation, the UV Mbh estimators are derived. Gon- 


sequently, a number of studies have been carried out to 
justify the use of Mg II or G IV-line based Mbh estima¬ 
tors comparing the masses from UV estimators to those 
from the optical. While some studies suggest that UV- 
line Mbh estimators are reasonably accurate, especially 
for Mg II, other studies point out a large scatter be¬ 
tween G IV-based measurements versus H/3-based mea¬ 
surements which can make the G IV -based Mbh values 
uncertain by a fa ctor of a few (e.g., iNetzer et al.ll2007l : 
iShen fc LuJl2012l hereafter SI2). It has b een noted that 
non- virialized motion (e.g., iDennevI 1201211 or extinction 
(e.g., lAssef et all 1201 111 , could severely modify the G IV 
line profile such that the Mbh cannot be reliably mea¬ 
sured from a simple virial equation, although it is con¬ 
troversial on the exact origin of the discrepant UV-based 
Mbh with that of optical. 

Purthermore, for the Mbh estimates to be valid, one 
also needs to justify the application of the low redshift 
Mbh estimators to high redshift, luminous quasars. Al¬ 
though one can expect that the Mbh estimators should 
not evolve in time in any significant way based on physi¬ 
cal ingredients of AGN models, this has not been tested 
at the high redshift, high luminosity regime. An ul¬ 
timate test would be to perform a reverberation map¬ 
ping study of high luminosity quasars at high redshift, 
but such a study would take decades to complete, since 
the variability timescale is long for luminous quasars 
and the cosmological time dilation makes it even longer. 
Another way, albeit less direct than the reverberation 
mapping method, would be to investigate the correla¬ 
tion between the line and continuum luminosities. At 
low redshift, the Ha or H/3 line luminosities are known 
to tightly correlate with the optical continuum luminos¬ 
ity (e.g., iGreene fc Hoi 1200^ . As the radiation energy 
L from an accretion disk increases, the distance to the 
broad line region increases as R oc ® from a sim¬ 
ple photo-ionization argument, or the energy flux of the 
radiation incident upon BLR clouds would be indepen¬ 
dent of the luminosity of the central p ower source, which 
has been confirined observationally (jKaspi et al.l 120001 
iKaspj et al.ll2005l : iBentz et al.ll20I,‘tl . A modification in 
the i?BLR”A relation at high redshift or high luminosity 
for example, would result in an increase or a decrease in 
the incident energy flux upon the BLR, thus the correla¬ 
tion between the line and continuum luminosities is likely 
to be modified. Interestingly, several studies suggest that 
the i?BLR — L relation is not valid for lu minous AGNs 
having massive black ho les with low spins (jLaor fc DavisI 
l20lHrWang et al.ll2014il . due to a decrease in the ionizing 
flux Lion caused by a decrease in the radiation tempera¬ 
ture of the accretion disk. Since the line luminosity Lnne 
is proportional to Lion, such theoretical expectations can 
be tested by examining the Liine - L 5100 relation at lu¬ 
minous end. 

In order to estimate the Mbh using an optical mass 
estimator and to test the universality of the scaling re¬ 
lations in the rest-optical for quasars at z > 3.5, spectro¬ 
scopic observation is necessary at A > 2.5 fim. This how¬ 
ever, is a very challenging task from the ground due to 
high thermal background at A > 2.5 /rm and atmospheric 
absorptions, limitin g such efforts to the study of Hff 
line at z<3.5 (e.g., Shemmer et aT]l2004 INetzer et ahl 
120071 : lAssef et S 1201 III . Recently, AKARI spectro¬ 
scopic observations have provided a breakthrough for the 
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study of the rest-optical spectra of distant objects, where 
its unique 2.5-5.0p,m coverage enables t he redshifted 
Hq lin e to be probed from 2 : = 3 to 6.5. lOvabu et al.l 
(j2007fl reported the AKARI detection of the redshifted 
Hq line of a q u asar R X J1759.4-1-6638 at 2 = 4.3, and 
iSedgwick et al.1 (j2013[ l investigated the star formation 
rate of submm galaxies at 2 > 3.5 based on Hq. 

With one of the AKARI mission programs (guaran¬ 
teed time) and also through several small open time pro¬ 
grams, we performed a spectroscopic study of quasars 
at high and low redshifts, with the aim to obtain the 
rest-optical spectra of high redshift quasars or the rest- 
frame 2.5-5.0^m spectra of low redshift AGNs. We call 
all these programs QSONG (Quasar Spectroscopic Ob¬ 
servations with NIR Grism) after the name of the mis¬ 
sion program, and here, we present the rest-optical spec¬ 
tral properties of 155 type-1 quasars from QSONG, along 
with scaling relations and Mbh estimates based on these 
spectra. 

The contents of the paper are constructed as follows. 
First, we use the spectro-photometric data including 
AKARI for 2 > 3 AGNs (section 2), and derive the 
continuum and line emission properties (section 3) in or¬ 
der to calibrate the Hq Mbh for the usage in high red¬ 
shift. We check the validity of the continuum and broad 
line luminosity relations and line width relations with 
respect to the local, update the mass equations for Ha, 
H/3, Mg II, and C IV, and compare the UV-optical Mbh’s 
(section 4). Finally, we discuss on the reliability of single 
epoch Mbh estimators at high redshift, and investigate 
the massive end evolution of Mbh in distant AGNs (sec¬ 
tion 5). Throughout this paper we adopt a flat ACDM 
cosmology with parameters of Hq = 70 kms“^ Mpc“^, 
Urn = 0.3, and Ha = 0.7. For the virial factor in the 
Mbh estimator, we adopt / = 5.1 ± 1.3 based on the 
calibration of the Mbh^ct relation using th e combined 
sampl e of AGNs and quiescent galaxies by IWoo et ^ 

(l2Ma . 


2. DATA 
2.1. Sample 


The majority of the data comes from the AKARI 
program QSONG, which is a two-year warm (phase- 
3) mission program consisting of ^ 900 approved point¬ 
ings, or ~ 150 hours of observations. The program is 
aimed to obtain the rest-optical spectra of high red¬ 
shift AGNs (I.Tun et al.l [MT 2 II or 2.5-5.0/xm spectra of 
low redshift AGNs (iKim et aH 1201511 containing Brack¬ 
ett and 3.3/xm Polycyclic aromatic hydrocarbon lines. 
The sample for this study is limited to the high red¬ 
shift, occupying 65% of the entire QSONG data. It 
is composed of optically luminous and spectroscopically 
confirmed typ e-1 quasars at -2 3, mostly out of SDSS 

(DR5 catalog . [Schneider et al.ll200^ and add i tional dis- 
fromlFan et al. 20001: lFa,n et al.ll2001l: iFan et akl 


coveries 


12001 IFan et al.l 12004 IFan et al.ll20n6irand APM-UKST 


( Storrie-Lombardi et al.l 119961 : iStorrie-Lombardi et al.1 
1200111 surv eys. Additional targets are from various 
references (iWebb et al.l 119881: iGregorv fc CondonI 119911: 
Henry et alT 1994 iGri ffith et al.l 1995t iKennefic k et al.1 


19951: Diorgovski et al.l Il998t iRengstorf et al.l 12004 

McGreer et al.l 1200611 . 


The targets are type-1 AGNs so that they allow the 



z 


Fig. 1.— 5100 A luminosity-redshift distribution of our AKARI 
observed quasars (black dots). The subsample with S/Nhq:> 2 
are pointed in red, and objects with rest-UV spectral coverage 
including the C IV emission are marked as blue squares. 

Mbh estimation from the broad line kinematics, through 
Hq/H/ 3 appearing within the AKARI near-IR (NIR) 
spectral coverage. In order to provide a minimal sensitiv¬ 
ity limit to the sample, we first considered the aperture 
size (68 cm) of the telescope and the restricted exposure 
time available at each sky position. After simulating the 
rest-optical spectra under the expected AKARI signal- 
to-noise ratio (S/N), the targets were chosen with 2 -band 
flux limits of ~18.5 and 19 AB magnitudes for bright 
and faint subsamples, respectively, with longer exposure 
time assigned to fainter targets for a clear line detection. 
Moreover, the targets were bounded in 3.3 < 2 < 6.4 so 
that the Ha emission and the surrounding continuum are 
placed within the 2.5-5.0/xm window of AKARI NIR 
spectroscopy. Without further constraints the targets 
were randomly selected in coordinates, redshift, and lu¬ 
minosity. We plot the distribution of ^-AsiqoA ™ Figure 
1 . 

Following the Hq observations of distant AGNs 
from AKARI NI R spectroscopy (|Ovabu et al.l l2007t 
lOvabu et al.ll200^ . our initial Hq survey of 14 quasars 
at 2 ~ 6 under th e Helium- cooled (phase-2) open time 
program HZQSO (Ilmll2010ll . demonstrated the feasibil¬ 
ity of AKARI observations in detecting the redshifted 
Hq emission. QSONG is essentially a phase-3 extension 
of the survey, for the purpose of vastly expanding the 
number of targets at the expense of warm phase sensi¬ 
tivity. Thus it probes a lower redshift distribution of 
quasars than HZQSO, with a peak at 2~4. In addi¬ 
tion, two more phase-3 open time programs, HQS02 
and DPQSO, were carried out either to push the red¬ 
shift limit of QSONG or to detect fainter optical lines 
(H/3 and [O III]), from deep exposures. Unfortunately, 
the Helium-dry observations led to significantly higher 
noise levels than expected, restricting the distinct scien¬ 
tific goals of the phase-3 programs that required better 
sensitivity. Therefore, we decided to merge all open time 
programs listed above under the scope of QSONG. 


2.2. Data acquisition 

We mostly used the NIR grism (NG) mode of the In¬ 


frared Gamera (IRC, lOnaka et al.U200H: lOhva ma et al 


l2007f) onboard the AKARI satellite ( Murakami et al 
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Fig. 2.— A sequential visualization of the data treatment additional to the pipeline processing. From the pipeline processed spectrum 
(left, ID at top, and 2D at bottom), hot pixels were rejected (center) before the spectral extraction within the aperture mask (red lines). 
The multiple pointings extracted (right, gray lines) were stacked (thick black line) with sigma clipping (red line indicates clipped data). 


TABLE 1 

Summary of Observations 


Program name 

Phase 

Observed period 

Mode 

Number of targets 

Observed pointings 

Rejected pointings 

HZQSO 

2 

Nov 2006-Aug 2007 

NG 

7 

26 

2 




NP 

5 

16 

7 

QSONG, HQS02, DPQSO 

3 

Jun 2008-Jan 2010 

NG 

147 

622 

74 




NP 

6 

11 

3 


Note. — NG and NP stand for NIR grism and prism modes, while a pointing is about 10 minutes long. For the QSONG program 
data, only the high redshift subsample is noted. The total number of targets is 165 excluding 22 rejected sources from confusion, 
problems in the spectra or data reduction (section 2.3). 10 objects are either duplicated or different in observing mode or program, 
yielding an effective total of 155 independent objects. 


I2007D . It offers a low, wavelength dependent spectral res¬ 
olution (i?), where i?=120 at 3.6 ^m. This corresponds 
to a velocity resolution of 2500 km s“^ in full width at 
half maximum (FWHM), sampled by a pixel scale of 
0.0097/rm in wavelength. The targets were placed in 
a 1' X 1' rectangular slit aperture to reduce source confu¬ 
sion. The wavelength dependence of R can be expressed 
as R = 12 0 (A/3.6 rtm) sinc e the dispersion is nearly a 
constant (ISakon et al.ll201§l . Meanwhile, a limited num¬ 
ber of NIR prism (NP, i?=19 at 3.5 p,m) observations 
were performed to better catch the fainter continuum 
and line luminosities. The angular pixel scale is 1.5" such 
that all targets are point-like in our probed redshifts. 

The observations were performed under the Astronom¬ 
ical Observation Template (AOT) mode of AOT04, typ¬ 
ical for spectroscopic observations. The number of NG 
pointings per target was normally 3-5 for QSONG, where 
one pointing observation corresponds to 355 or 400 sec 
on-source exposure. The number of pointings were de¬ 
termined based on the 2 -band flux, generally set to be 
smaller for phase-2 and larger for phase-3 open time 
programs. The NP observations were performed with 
usually 1-2 pointings. By the time of termination of the 
satellite mission, the QSONG program was 85% com¬ 
plete with 144 high redshift quasars observed. The open 
time programs were complete before the satellite lifetime, 
adding another 33 targets. Table 1 summarizes the ob¬ 
servations. 

To supplement the NIR spectra, we compiled the op¬ 
tical spectra of the AKARI quasars from the Sloan 
Digital Sky Survey (SDSS) database (DRIO including 
both the SDSS-I /SDSS-II and the SDSS-III BOSS data, 
lAhn et ^ |2014[) , and from observ ations of APM-UKST 
quasa rs (IStorrie-Lombardi et al.l [ 1996: iPeroux et HI 
[20?n^ and 00000-26 (ISchneider et all Il989ir in order 


to estimate the C IV line based Mbh- Also, we col¬ 
lected broad-band photometric data from optical to 
mid-infrared (MIR) imaging, for the calculation of the 
rest-frame UV-optical continuum luminosity of AGNs 
through SED fitting (section 3.2). The data includes 
SDSS DR9, 2MA SS PSC, UKIDSS DRIP, and WISE 
AllWISE releases (lAhn et al.ll2012 USkrutskie et al.l[2006t 
[Lawrence et al.l 120071 : IWright et al.l 1201011 . and ex- 


isting Pan-S TAR RS, Spitzer , and AKARI i magin; 
llHines et al. I200d Jiang et al.l 120061: lOvabu et al.ll200L . 
iJiang et al.l 20101 Leioski et al.l l20lT The Galactic 
extinction is corrected for these spectro-photometric 
data, assuming the total-to-selective ex tinction ratio of 
Rv=3.1 and using the corrected form (IBonifacio et al 

a to the extinction map values of iSchlegel et al 
. The photometric measurement schemes are dif¬ 
ferent in each survey data, such that host galaxy contam¬ 
ination may not be well subtracted. We keep the diverse 
magnitude type and imaging resolution of each survey 
data however, as the high luminosity AGNs yield a com¬ 
patible set of magnitudes dom inated by the central AGN 
contribution (|Jun fc Imll2013ll . The optical spectra and 
the multi-wavelength imaging data are outlined in Table 


2 . 


2.3. Data reduction 

The data were reduced using the automated IDL 
pipeline package IRC_SPE CRED (versions 20 110114, 
20111121 for phase-2 and 3. IOhvama et al.ll2007[l . where 
pre-processing (dark, linearity, flat corrections), image 
registration and coaddition, flux and wavelength cali¬ 
bration, astrometry, spectral extraction, and aperture 
correction were the main tasks involved. The standard 
pipeline configuration was adopted, except for the us¬ 
age of short exposure when the image taken for regis- 
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tering subframes of spectral data were contaminated by 
saturated stars. This procedure considerably improved 
the registration of 2-D spectra in both spatial and wave¬ 
length directions. In addition, astrometry of the refer¬ 
ence image was upgraded using the 2MASS point source 
catalog to better extract the faint NIR spectra. Indeed, 
the zeroth order positions of the spectra and the 2MASS 
coordinates were visually well aligned for spectral extrac¬ 
tion. 

Due to the increased number of hot pixels and the 
background level in the phase-3 data, a non-negligible 
number of bad pixels remained in the reduced spectral 
data even after applying the IRC_ SPECKED pipeline. 
To remove the remaining hot pixels, a further data re¬ 
duction step was taken to obtain a cleaner set of spec¬ 
tra, as depicted in Figure 2. After running the spectro¬ 
scopic pipeline we subtracted the rem aining hot pixels 
using L.A.Cosmic (Ivan DokkumI [20Mll with a threshold 
of 2.5(7, and combined the 1-D spectra with a 2.5-sigma 
clipping. This threshold was chosen by visually inspect¬ 
ing the reduced spectra so that the chosen threshold re¬ 
moves the spiky hot pixels efficiently without affecting 
emission lines. Because the typical FWHM of the broad 
lines in our AGNs is broader than the spectral resolution 
of NG (2500kms“^) but not that of NP (15,800kms“^ 
at 3.5/xm), the hot pixel rejection and the sigma clip¬ 
ping were applied on the NG data only, while the under¬ 
sampled NP emission spectra were kept unchanged. We 
clipped 2.7% and 2.9% of the NG spectral data through 
the hot pixel reduction procedure and the combining pro¬ 
cess respectively. 

We extracted the 1-D spectra of 3 pixel width for 
the phase-2 data to maximize S/N, but the width was 
widened to 5 pixels for the phase-3 data since the noisier 
spectra made it difficult to determine the center of the 
object spectrum. Host galaxy contamination is negligi¬ 
ble at the bright luminosities of the sample quasars (Sll), 
enabling flexible extraction widths. Aperture corrections 
were automatically carried out from the pipeline to de¬ 
rive the total flux, for given respective extraction widths. 
The pipeline did a fair job of placing the extraction aper¬ 
ture on the right location, but visual inspection showed 
it necessary to make a —1 pixel shift in spatial direction 
for 90% of the sample. For 10% of the cases, a different 
shift of —2 to 1 pixels was necessary. The wavelength ze- 
ropoints were determined from the pipeline taking into 
account the satellite attitude drift and sub-pixel coordi¬ 
nate rounding effects, and we did not apply any further 
correction as the zeropoints were confined within a 0.5 
pixel scatter. 

Multiple pointings of the 1-D spectra were stacked for 
each object, where the NG spectra were interpolated to 
a fixed wavelength grid, flux averaged, and error rescaled 
assuming Poisson error statistics. Multiple pointings of 
NP spectra were stacked without modifying the individ¬ 
ual spectrum, due to their poor resolution. We pro¬ 
vided secondary flux calibration to the stacked spectra 
by integrating the AKARI fluxes over the WISE filter 
response curves, to match the W1 and W2 fluxes to¬ 
gether by a constant additive correction. The average 
and rms scatter of the corrections are — 0.01 ± 0.19mJy. 
Out of 675 pointing observations, we used 589 pointings 
from 165 objects, since some of the pointing observations 
were not usable due to contamination of the object spec- 


TABLE 2 

Supplementary Data 


Name 

Wavelength 

N 

Exposure 

Reference 

Spectra 

SDSS 

3800-9200 A 

111 

> 45m 

1 

BOSS 

3650-10400 A 

98 

> 45m 

1 

APM-UKST 

~3500-9000 A 

16 

15-60m 

2,3 

Hale 

4500-9000 A 

1 

30m 

4 

Photometry 

SDSS 

ugriz 

134 

54s 

5 

Pan-STARRS 

zy 

14 

6m 

6 

2MASS 

JHK 

76 

8s 

7 

UKIDSS 

YJHK 

34 

40s 

8 

WISE 

W1-W4 

161 

200s 

9 

Spitzer 

IRAC, MIPS 24fj.m 

27 

17-23m 

6,10,11,12 

AKARI 

IRC 2-11/im 

1 

- 

13 


Note. — N is the number of matches to the 165 
AKARI objects in Table 1. Exposure t imes are typical 
val ues. The references number ed ar e 1. lAhn et aIjl2Q1 4l: 
2. ISto rrie-Lornbardi et al 1119961: 3. IPeroux et al.l I20Q1I: 4. 
Schneid er et al.| 119891; 5. lAhn et al.l 120121: 6. ILeipski ^ al.l 

20141: 7. ISkrutskie et al.l 1200 61: 8. ILawrence et al.l I20Q7I: 9 . 

Wright e t al.ll20 10l: 10 . IHines et al.ll2006| ; 11. IJiang et al.ll2Q06l : 
12. IJian^^^aLil2010l : 13. IOvabi^^^.1120071 . 



Fig. 3.— The S/Nhq- distribution of the 160 AKARI observed 
quasars with available S/N (section 3.1). The S/Nhck >2 and >3 
cuts are marked in dotted and dashed lines, which are given to 
limit the measurement of Lna and FWHMHa respectively (Figure 
5). 



Fig. 4.— Examples of the rest-frame Ho emission fitting. On 
top of the spectra (thin line), the resolution (3 pixel) matched 
data (dots) and errors, the best-fit to the continuum and Ha line 
emission (thick line), and FWHM are indicated. On the figures, the 
observation mode, the name of object, the redshift and S/N of the 
Ha emission, the number of AKARI pointings, Npt, Ha luminosity 
(ergss“^) and FWHM (kms“^) are printed. The (a) NG and (b) 
NP observations performed for BR J0006-6208 gives an idea of the 
enhanced resolution and sensitivity of each spectroscopic mode. 
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Fig. 5.— Top rows: An example fit to the resolution degraded SDSS spectrum to test the reliability of Mbh from AKARI observations 
(left), where the continuum subtracted original SDSS spectrum (gray), and the AKARI resolution matched spectrum (black dots) with 
fit to the simulated data (thick line) are plotted. The broad FWHMjja values of 2 < 0.3, T 5100 > 10“^® ergs s“^ SDSS quasars from Sll 
are compared to our fitting result of the simulated spectrum (center). Offsets from a 1-1 relation and intrinsic scatter are displayed. For 
the comparison of FWHMs (center), we divide the sample into FWHMjja larger or smaller than 2500 kms“^ (filled and open dots) and 
only use the FWHMho > 2500 kms“'^ data. The red highlighted symbol represents the example on the leftmost panel. Likewise, the Ho 
luminosity from Sll and our single Gaussian fit to the simulated spectra, are compared (right). Bottom rows: we also test the effect of 
low S/N to the fitting by adding a set of Gaussian random noise to the resolution degraded spectrum (left). The ratio of our FWHMjjq, 
Thu measurements to that from Sll are plotted along the S/Nho (center and right). For the comparison of FWHMs (center), we remove 
49 data points from the plot with AFWHMjja = 0 usually at FWHMjjq < 2500 kms~^. The mean and 1 a offsets are shown in black and 
gray lines. When calculating the mean and intrinsic scatter of the quantities, only the S/Nhci >3, FWHMhq >2500 kms“^ data are used 
for the comparison of FWHMjjq, and the S/Njja > 2 and any FWHMjjo data for the Lua- 


trum by adjacent sources. We excluded frames from the 
analysis when there was a source that is brighter than 
the target and its distance from the target is less than 
the FWHM in spatial direction. In rare occasions, the 
reduction pipeline did not run properly and such data 
were not used. In Figure 3, we plot the histogram of 
S/Nhci, the S/N of the Ha emission line measured within 
± FWHMho from the line center. 

3. ANALYSIS 
3.1. Spectral fitting 

We modeled the rest-frame 5500-7500 A spectra of 
our AGNs, as a sum of the power-law continuum 
/a = (where f,j oc and the Gaussian Ha 

emission components. We did not attempt to fit the 
relatively weak emission features far (He I, [O I]), or 
near the Ha ([N II] and [S II] doublets, Fe II complex), 
as they were not detectable under the AKARI spectro¬ 
scopic resolution and sensitivity . Likewise, the H/3 and 
[O III] lines were too weak to be detected in most cases 
and were not fitted accordingly. Examples of the fit¬ 
ting are given in Figure 4, with fitted parameters shown 
on each panel. We found that 23% of the sample show 
spiky emission/absorption features around the Ha line 
from low S/N spectra, which were manually masked out. 
The Ha line was modeled as a single broad Gaussian 
with observed FWHMs of 2500-10,000km s“^, for the 


inability to clearly resolve the narrow or multiple broad 
components at 120. The Ha line center was set free 
within ±2500km s“^ to the UV line-based redshift from 
references in Table 2. However, we found four (SDSS 
J143835.95±431459.2, SDSS J142243.02±441721.4, HR 
J0307-4945, and SDSS J150654.54±522004.6) exceptions 
whose Ha line centers were significantly redshifted from 
the UV-line based redshifts. In these cases, the veloc¬ 
ity range was loosened to ± 10,000 km s“^, where the Ha 
showed velocity shifts of 3600-5600km s“^. Next, the 
measured broad line width FWHMobs, was subtracted 
by the instrumental resolution FWHMjns (section 2.2) in 
quadrature, to obtain the intrinsic line width FWHM = 
A/(FWHMobs)2 - (FWHMi„,)± 

Extracting the broad emission line luminosity and 
width is important in accurately estimating the Mbh 
of AGNs. The limited AKARI sensitivity (Figure 3) 
and spectral resolution could produce systematic bias in 
the measurement of line p arameters such as th e line lu¬ 
minosity and width (e.g., iDennev et al.l[2?)0^ . There¬ 
fore, we investigated how the low resolution, low S/N 
spectra systematically affect the results of the spectral 
fitting. This was done by running Monte Carlo simu- 
lations on a set of lum inous SDSS quasar spectra (DR7, 

I Schneider et al.ll2010D of i? ~ 2000 to mimic the quality of 
AKARI spectra. Given that the host galaxy contamina¬ 
tion to the quasar spectrum is negligible at L 5100 > 10"^® 





























7 



0.55 0.60 0.65 0.70 0.75 0.55 0.60 0.65 0.70 0.75 0.55 0.60 0.65 0.70 0.75 0.55 0.60 0.65 0.70 0.75 



0.55 0.60 0.65 0.70 0.75 0.55 0.60 0.65 0.70 0.75 0.55 0.60 0.65 0.70 0.75 0.55 0.60 0.65 0.70 0.75 



0.55 0.60 0.65 0.70 0.75 0.55 0.60 0.65 0.70 0.75 0.55 0.60 0.65 0.70 0.75 0.55 0.60 0.65 0.70 0.75 



Fig. 6 .— Rest-frame Ho; spectral fitting of selected objects with S/Nhq; >3, sorted by the Ha redshift. The data point symbols and 
colors follow the meaning of Figure 4. When there were contaminations to the Ha emission, we masked out the region (red). 


ergss ^ (Shen et al. 2011, hereafter Sll), we collected 
15, isioo > 10^^ ergss”^, type-1 quasar spectra with con¬ 
tinuum sensitivity of S/N>20, where the Ha emission 
region is present (z < 0.3) and well fit {xt < 2) from Sll. 

We smoothed the SDSS spectra with a Gaussian func¬ 
tion to match the AKARI NG resolution, i?=120, and 
rebinned the data to match the 3 pixel per resolution 
sampling of the AKARI spectra. First, we looked into 
the question of how the results get affected with a single 
Gaussian fit to the emission line in the low resolution of 
the AKARI data. For this, without adding extra noise, 
we followed the same method to measure the Ha line 
FWHM and luminosity as for the AKARI spectra (e.g.. 
Figure 5a). In Figures 5b and 5c we compare the fit¬ 
ted parameters FWHMna and Lna from the smoothed 


and binned spectra, to that of the measurement from 
SII. We find that the FWHM at >2500kms“^ and 
the line luminosity from the degraded resolution spec¬ 
tra, are remarkably consistent with Sll within ^ 0.04 dex 
offset and intrinsic scatter Uint, where — 

f{xi)Y - Ayf - l3‘^Ax‘f]/{N - 1) for f{x) = a + Px 
and a set of N data points {xi, yi) with measurement er¬ 
rors (Axi, Aj/i). This could bias the Mbh measurements 
up to ^0.1 dex in offset and a-mt when following the 
Mbh ^ x FWHM^ behavior, but this is smaller than 
the typical Mbh measurement uncertainty (e.g., Figure 
14c). Exceptions to the consistency between the simu¬ 
lated and observed parameters are the line width mea¬ 
surements at FWHMhq; <2500 kms“^, where the sim¬ 
ulated FWHM values fall below the extrapolated linear 
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SDSS J083700.82 + 350550.2 
z = 3.31. S/N = 25.2 
logL,350=46.92. FWHMciv=4044 


x!=1-4 



SDSS J142243.02+441721.4 
z = 3.55, S/N = 23.1 
logL,350=47.04, FWHMc,v= 12353 


SDSS J144144.76 + 472003.2 
z = 3.63. S/N = 8.5 
logL,350=46.71, FWHMciv=3313 
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Fig. 7.— Rest-frame C IV spectral fitting of the identical objects in Figure 6 plotted at the same relative location, following the format 
of the plotted data and printed numbers. Additionally shown are the double Gaussian fit (gray) to the C IV above the continuum, and 
the 1350 A monochromatic flux marked on its wavelength (blue dots), with extrapolated values indicated (dashed lines). The spectra are 
shown when they are available from the listed references in Table 2. When there were contaminations around the C IV emission, we masked 
out the region (blue). The C IV spectra displayed are smoothed down to i? = 500, to highlight the spectral features better. 


relation. Therefore, we give a FWHMhq: > 2500kms“^ 
limit to our AKARI data to restrict the sample with less 
biased line width measurements. Meanwhile, the narrow 
Ha and [N II] doublet luminosities of the fiducial SDSS 
AGNs add up to the broad Lho, by less than 0.01 dex 
altogether. The weakness of the narrow emission lines in 
luminous AGNs guarantees that the narrow line contam¬ 
inations to the degraded resolution spectra are negligible. 

Second, to investigate the effect of low S/N to the fit¬ 
ted results, we added a set of random Gaussian noises 
on top of the degraded resolution spectra. Each SDSS 
spectrum was repeatedly simulated 30 times adding the 
random errors, for the S/Nhq, to be distributed down to 


the level of AKARI S/N. Again, we followed the same 
method to measure the Ha line FWHM (only for the 
FWHMho > 2500km s“^ objects from Figure 5b) and lu¬ 
minosity as for the AKARI spectra (e.g., Figure 5d). In 
Figures 5e and 5f we plot the ratios of FWHMho and Lho, 
from the degraded resolution and S/N spectra, to that 
from Sll, along S/Nhq- Overall, we find both FWHMho 
and Lhq to be within 0.02 dex in the mean offset to the 
values before adding the noises (Figure 5b and 5c), at 
FWHMna >2500kms-7 S/NHa>3 for the FWHMh^ 
and S/Nho > 2 for the Lho;- Also, the (Tint in FWHMho 
and Lho under low S/N are within 0.03 dex to those of 
the noise-free, effectively unchanging the bias in the Mbh 
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SDSS JO837OO.02 + 35O55O.2 
z = 3.31, a = -O.OB, logLgjoo=46.62 

^oS^i350“46.94 







^re.l (fim) 





^r.si (Mm) >'r„i (Mm) X„,, (;j.m) 

Fig. 8 .— Rest-frame UV-NIR broad-band SED of the objects in Figure 6 plotted at the same relative location. The figure shows the 
observed data points (filled circles) and WISE 2 a upper limits (arrows), model fits of the accretion disk (blue line) and the T=1250K 
dust components (red curve). Also, the monochromatic 1350, 3000, and 5100 A fluxes are drawn (open circles) when available. 


at a similar ^O.ldex. We note that below the S/N or 
FWHM limit the fitted quantities systematically diverge 
from those calculated with the original spectrum. There¬ 
fore, we conclude that for a given selection of data nei¬ 
ther the poor resolution nor sensitivity biases the fitted 
results by greater than ^O.ldex level of systematic off¬ 
set or scatter in Mbh, and we give the corresponding 
FWHM and S/N cuts to the AKARI data. 

Having tested the reliability of fitted quantities under 
possible systematic biases, we come back to the fitting 
of the AKARI data and find the fit to converge for 160 
out of 165 objects. Five failures show noisy spectra near 
the Ha and were removed. The goodness of the spec¬ 
tral fitting is quantified as the reduced chi-square, and 
it reaches down to = 0-33 in median. The fraction of 


AGNs passing the reliability limit for the FWHM mea¬ 
surements (S/Nho > 3) are 67% (N=8) and 25 % (N=37) 
for the phase-2 and 3 data respectively. Also, we flagged 
the five NG objects of FWHMjnt < 2500 km s“^ with an 
upper limit of 2500km s“^, and put upper limits on the 
measured FWHM of the four NP sources (e.g., Figure 
4b). In total, there are 43 FWHMhq measurements in¬ 
cluding seven out of nine upper limits, where the ex¬ 
cluded two upper limits are NP measurements with over¬ 
lapping coverage in NG. 

Next, we computed the Ha line luminosity Lua and 
the 5100 A continuum luminosity T 5100 from the spectra, 
converting the measured rest-frame fluxes using the lu¬ 
minosity distance and assuming isotropic radiation. The 
Ljja was derived from the Gaussian fit to the observed 
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Fig. 9. — Comparison of spectroscopically and photometrically derived monochromatic luminosities, for L 1350 (left) and L 5100 (right)- 
The AKARI NG/NP data points for the I /5100 are highlighted with red squares and yellow triangles, respectively. Offsets to the luminosities 
from a 1-1 relation is displayed. 


flux. The reliability limit for Lhq (S/Nhq > 2) is satisfied 
for 45% (N=72) of the data, while the S/Nhq < 2 spec¬ 
tra were provided with 2 a upper limits from their given 
noise levels collected within ± 4000 km s“^ of the Hct line 
center. Meanwhile, the T 5100 was calculated from the av¬ 
erage of the rest-frame fluxes at 5000-5200 A to reduce 
the measurement uncertainty. The T 5100 measurements 
were kept only when S/N 5 ioo> 2 , and the rest of the 
data were given with 2 a upper limits alike Lho- We 
have less T 5100 measurements than Ljja where only 25 % 
(N=41) meet S/N5ioo>2, not to mention the limited 
number of spectra (53%, N=88) covering the rest-frame 
5100 A. Thus, we also derived the T 5100 alternatively by 
the photometric SED fitting (section 3.2). 

In addition, we fitted the C IV region (rest-frame 1445- 
1705 A) of 121 objects with a SED model containing a 
power law component and double broad Gaussians to 
model the C IV emission. Also, a single broad Gau ssian 
was used to fit the 1600 A feature (jLaor et al.lll994D . and 
the He II and O III] around 1650 A altogether since these 
emission are blended but relatively detached from the 
G IV. This component was not regarded as a part of the 
G IV, co nsistent with the p revious studies (e.g,. prescrip¬ 
tion A of lAssef et al.ll201ll : S12). We do not subtract the 
broad Fe II complex around the G IV emission, as it does 
not change the FWHMciv meaningfully (Sll). The op¬ 
tical spectra were fitted after carefully masking out the 
absorption features around the C IV line for 25 % of the 
spectra. Meanwhile, 16 spectra without error informa¬ 
tion were fitted assuming the flux error is uniform at all 
wavelengths, and the rms scatter of the best-fit solution 
is chosen to be the flux error afterward. Out of five spec¬ 
tra with severe broad absorption line (BAL) features, we 
fitted the G III]A1908 and used its line width as an ef¬ 
fective FWHMciv surrogate (S12) for two objects, while 
excluding the remaining three objects from the UV line 
analysis. In total, we derived FWHMcjv and T 1350 for 
118 objects. The T 1350 and its error were calculated from 
the average of the rest-frame 1350 ± 15 A fluxes to avoid 
contamination from narrow absorption, while for seven 
BAL quasars we extrapolated the continuum around the 
G IV emission to 1350 A. When spectra of an object 
were available from both SDSS-I/SDSS-II and BOSS, we 
performed the fit to the spectra from each dataset sepa¬ 


rately, and took the average of the parameter values from 
the independent fits. We plot examples of the spectral 
fitting of the Ha region in Figure 6, and of the C IV in 
Figure 7. 

3.2. Broad-band SED fitting 

The photometry datasets in Table 2 cover a wide wave¬ 
length range in broad-band filters from u-band through 
24/rm, thus we fitted the broad-band SEDs to provide 
further information on the AGN continuum luminosities. 
Under a photometric sensitivity limit of S/N > 5 for the 
observed optical-NIR and S/N > 2 in the MIR, and fur¬ 
ther rejecting the 2MASS data with a single filter de¬ 
tection, we modeled the SED in the rest-frame 0.3-5 ruy 
including 5 data points on average, as a sum of a power 
law continuum an d a black body e mission from hot dust 
ofT=1250K (e.g.. l,lun fc Iml[2Cil^ . For 53% of the case 
(N=87) for which there were no rest-frame 0.3-0.6/im 
data points available, we used the average continuum 
slope a = -0.08 of luminous SDSS quasars (lJun fc Iml 
120131) . The uncertainty in fixing the continuum slope 
was tested from objects that cover the rest-frame 0.3- 
0.6/rm, and comparing the L 5100 with and without fix¬ 
ing the a value. The test yields the L 5100 to be offset 
by — 0.01 ± 0.04 dex when a is fixed, compared to when 
a is free. Since the offsets are small, we find our method 
to fix the a when missing the photometric coverage near 
the rest-5100A, to be reliable in tracing the Lsioo- In 
addition, considering the filt er bandwidths and the AGN 
line equivalent widths from IVanden Berk et ^ (|200lD , 
we find the Ha to be the only line that meaningfully 
contributes to the rest-optical photometry over the con¬ 
tinuum emission (by >0.03dex). Thus, we removed the 
data point enclosing the Ha emission when the con¬ 
taining that data point became larger than that without. 
Through this procedure we obtained 164 photometrically 
derived L 510 O) while removing one object without any de¬ 
tections in the rest-frame 0.3-5 fj,m under our sensitivity 
limit. Examples of the broad-band SED fitting are shown 
in Figure 8. The reduced chi-square values have a me¬ 
dian of = 3.6, is acceptable given the simplification of 
the SED model that does not take into account the emis¬ 
sion line features, and the general agreement of the fit to 
the data demonstrated in Figure 8. 

Meanwhile, we interpolated the broad-band SED 
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TABLE 3 

CONTINUM AND LiNE BASED PROPERTIES OF THE SAMPLE 


Name 

-^ref 

^Ha 

log L 1350 

log ±5100 

log ±Ha 

FWHMs.Civ 

FWHM3,hc< 

log Mbh, Civ 

log Mbh,H a 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 


(7) 

(8) 

(9) 

(10) 

SDSS J000239.39-I-255034.8 

5.80 

5.79 

99.00 ±99.00 

46.56 ±0.02 

45.14 ± 

0.11 

99.00 ±99.00 

2.50±-1.00 

99.00 ±99.00 

9.25±-1.00 

Q 0000-26 

4.10 

4.11 

47.40 ±0.01 

47.13 ±0.01 

45.48 ± 

0.13 

4.33 ±0.28 

2.50±-1.00 

9.88 ±0.22 

9.56±-1.00 

SDSS J000552.34-000655.8 

5.85 

5.85 

99.00 ±99.00 

46.04 ±0.08 

44.75 ± 

0.32 

99.00 ±99.00 

99.00 ±99.00 

99.00 ±99.00 

99.00 ±99.00 

BR J0006-6208 

4.45 

4.49 

46.90 ±0.02 

46.71 ±0.00 

45.46 ± 

0.06 

11.33 ±1.33 

2.85 ±0.85 

10.48 ±0.24 

9.46 ±0.31 

SDSS J001115.23-1-144601.8 

4.97 

99.00 

99.00 ±99.00 

46.78 ±0.02 

47.97 ±- 

-1.00 

99.00 ±99.00 

99.00 ±99.00 

99.00 ±99.00 

99.00 ±99.00 


Note. — Catalog of the properties derived for the AKARI quasars, sorted by right ascension. Column 1: Target name; Column 2: Redshift 
from references; Column 3: Redshift measured from Hen; Column 4: 1350 A luminosity and its uncertainty; Column 5: 5100 A luminosity and its 
uncertainty; Column 6: Hen luminosity and its uncertainty; Column 7: FWHM of the C iv line and its uncertainty; Column 8: FWHM of the 
Ha line and its uncertainty; Column 9: Mbh from the C iv line and its uncertainty; Column 10: Mbh from the Ha line and its uncertainty. 
The units for L, FWHM, and Mbh ergss”^, 1000kms~^, and Mq. Columns 9 and 10 are from Equations (10) and (7), respectively. Empty 
parameters are entered as 99 and upper limits are given with errors of -1. 

(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form 
and content.) 
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Fig. 10. — Composite spectra from AKARI observations, (a) The spectra within 0.4-1/im (top), and the number of spectra used to 
construct the composite, Nc, plotted against wavelength (bottom). The mean Nc and the continuum slope from the 0.5-0.9/im region are 
printed, (b) Ha region fit of the composite spectra, (c) region fit of the composite spectra. The sum of the continuum, emission, 
and the Ee II complex (thick black line), and the continuum (gray line) are overplotted on the data, while the H/3/[0 Hi] and the Ee H 
emission components are separately plotted below the spectrum in black and gray lines. 


around the rest-frame 1350 and 3000 A to obtain L 1350 
and Laooo- For this, we used objects with more than 
two data points in the rest-frame 500-2500 A for L 1350 
or 2000-6000 A for Lsooo- The interpolation is done lin¬ 
early to the data points, and we obtain 137 T 1350 and 
47 T 3000 values. The continuum luminosities derived 
by photometric and spectroscopic methods roughly agree 
with each other as shown in Figure 9 for L 1350 and T 510 O) 
though there are not enough data points (N=2) to plot 
for ^ 3000 - Likewise to L 5100 , we calculated the level of 
C IV or Mg II line contamination to the L 1350 and L 3000 
from broad-band photometry. The C IV and Mg II el¬ 
evates the observed broad-band flux by up to 0.06 and 
0.04 dex. Thus, it is possible that the photometry em¬ 
bracing the broad UV emission is boosted by more than 
the typical measurement error in L 1350 and Lsoooj which 
are 0.02 and 0.03 dex, respectively. Between the spec¬ 


troscopically and photometrically derived continuum lu¬ 
minosities, we will use in the following discussion the 
spectroscopically derived L 1350 and the photometrically 
derived L 3000 and Lsioo- We do so because the line con¬ 
taminations near 1350 A through the broad-band pho¬ 
tometry can be significantly larger than the measurement 
uncertainty, while it is not so around the Mg II and Ha 
lines. The large error in L 5100 from AKARI spectra is 
also another reason why we opt to use L 5100 from the 
broad-band SED fitting. The Htted properties from this 
section, and the Mbh to appear in section 4, are listed 
in Table 3. For objects with both NG and NP observa¬ 
tions, we use the zho and FWHMhq, from the NG and 
the Lhq: from the NP when S/Nho > 3, while we list only 
the values from the NG otherwise. 

4. RESULTS 
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Fig. 11.— The Tsioo^l'Ha relation of AGNs (top), and its projection on the Lsioo^l^Ha/Tsioo (bottom), from combined references. 
The references abbreviated on the plot are summarized in Table 4, where all measurements from references are converted to our adopted 
cosmology. We limited the G05 data to Lsioo < ergss“^ to avoid overlap with the Sll data. The Tsioo^Tho relation from G05, S12, 

and J15 are shown in dotted, dashed, and solid lines respectively. The zeropoint offset and rms scatter of each literature data with respect 
to our relation are denoted as A^pt. The intrinsic scatter (cinti dex) of the entire data with respect to the G05, S12, and J15 relations 
(top), and of the data divided by L 5100 = 10^“^ ® ergss“^ (bottom), are shown. The T 5100 in M13 and J15 are from photometric SED fitting 
(section 3.2), while the rest are spectroscopically derived from each reference. The Lua are from broad emission for all references but for 
the G05, 307, and Sll data, where the narrow component is included. 


4.1. Composite spectra 

To investigate the overall rest-optical spectral proper¬ 
ties of the sample, we construct the composite AKARI 
spectra. Out of 154 objects observed with NG, 127 are 
used for the composite construction after removing 27 
spectra due to a mild level of confusion from neighbor 
source spectra, negative continuum levels, or strong fluc¬ 
tuations near the Ha due to warm pixels. The composite 
is constructed by normalizing the spectra at 5100 A and 
taking their error weighted mean to maximize the S/N. 
Each spectral flux and error were de-redshifted, and re¬ 
binned to a common wavelength scale of 18 A per bin 
which is equal to that of the AKARI at rest-frame Ha. 
In Figures lOa-lOc, we plot the composite spectrum, and 
zoomed-in fit to the Ha and H/3 regions, respectively. 
The Ha emission is prominent in the composite spec¬ 
trum, and the spiky feature at 8400 A is an artifact from 
a single spectrum with high S/N. 

We determine the continuum slope from the 0.5-0.9 /rm 
region from Figure 10a, where the number of spectra 
used to construct the composite exceeds 60. The slope 
a = —0. 52 ±0.06 (where /t, oc v °‘) is close to a = — (0.37- 
0.48) of iGlikman et all (|2006D determined through the 
composite of local luminous quasars at similar wave¬ 
lengths, indicating a similarity in the rest-optical con¬ 


tinuum shape of luminous type-1 quasars with respect 
to redshift. Interestingly, we detect a sign of the H/3 
emission fro m Figure lOc. The H/3 region was fit¬ 
ted with the iBoroson k. GreenI (119921) Fe II template, 
following the method of iShen et al.l (l2008f l. We find 
FHa/Fn/S = 4.5 ± 1.6, which is roughly consistent with 
3.6 ±1.4 from luminous z ~ 2 quasars (S12) or the model 
broad Balmer line decrement of AGNs at T= 10,000- 
12 ,000K, LHa/iH /3 = 3.6-8.8 (lOsterbrockllTMl) . Since 
the Balmer decrement value and the S/Nhq of each 
AKARI spectrum suggest that the strongest H/3 in our 
individual spectrum would appear as S/NH,a=l“2, we 
do not expect the H/3 emission to be individually de¬ 
tected, consistent with the visual inspection in section 

3.1. Apart from this, we do not find hint of other emis¬ 
sion lines in the composite spectra. 

4.2. Luminosity and Line Width Scaling Relations 

The derivation of continuum and line luminosities 
for distant, luminous quasars allows us to examine the 
universality of the luminosity scaling relations over a 
wide range of redshifts and luminosities. Starting from 
the ± 5 ioo“AHa relation, we plot in Figure 11 our de- 
rived data po i nts a nd the data tak en from literatures 


(IGreene & Hoi 1200,51 hereafter G05; IShansr et al.l 

2007 

hereafter S07; Sll; 

Ho et al.ll2012l:IMatsuoka et al.l 

2013 
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Fig. 12.— The Tsioo^-^isso I/ 5100 —L 3000 relations of AGNs from combined references. The references abbreviated on the plot are 

summarized in Table 4, where all measurements from references are converted to our adopted cosmology. The T 1350 of the N07 sample 
are searched from Sll. The T 5100 , T 3000 , and Z /1350 are derived from spectra, except for the L 5100 , ^3000 from J15 where they are from 
photometric SED fitting (section 3.2). We assign a modest 10% error for L 5100 of t-Ee NOT data and L 1350 , T 5100 of the All data, and 
20% error for the L 1350 , L 3000 , and L 5100 from D09 data, from visual inspection of their spectra. We removed the highly variable object 
3C 390.3 from the B13 data, and additional two objects in SOT data that overlap with B13. The I/5ioo~-^'i350? T5ioo~-^'3000 relations from 
S12 and J15 are shown in dotted and solid lines respectively. The zeropoint offset and rms scatter of each literature data with respect to 
our relations, are indicated. 


TABLE 4 

Dynamic Range of References 


Reference 

Z 

log L 5100 (ergss 1 ) 

G05 

<0.35 

41.7-45.0 

NOT 

2.3-3.5 

45.2-46.3 

SOT 

0.08-0.40 

44.4-45.8 

M08 

0.35-0.3T 

43.6-44.4 

D09 

1.1-2.2 

46.1-46.7 

W09 

0.00-0.16 

42.0-45.9 

All 

1.4-3.6 

44.8-46.5 

Sll 

0.08-0.80 

43.8-46.4 

H12 

1.4-1.5 

45.2-46.0 

S12 

1.5-2.2 

45.4-47.0 

B13 

0.00-0.29 

42.0-45.9 

M13 

O.T-l.T 

44.8-45.5 

P13 

0.01-0.23 

42.6-45.9 

J15 

3.3-6.2 

46.0-47.2 

Note. — The abbreviated references are IGreene &: Hoi 

20051 (G051: INetzer et al.l 120071 (NOT): 

Shane et al.ll2007| (^SOTh 

McGill et al.ll200S rMOSh iDietrich et al. 

120091 TDOOh IWane et al.l 

20081 fW09i; lAssef et 

al.l 120111 rAllh 

IShen et al.l 120111 fSlli: 

Ho et al.l 12013 fH121: 

Shen & Liul 120121 (S12): IBentz et al.l 120131 

fB13h IMatsuoka et al. 

120131 fM131: IPark et al.ll2013l ('P13'): and 

this work (J15). 



S12) that cover a 

range of A 5100 and 2 , as summa 

rized in Table 4. 

Our AKARI data extends the re 

lation at the high redshift (z > 3.3) and high luminos 


ity end (L 5100 > 10“^® ergss ^). To minimize the host 
galaxy contribution to the AGN luminosities, we chose 
AGNs with host contamination <20% in L 5100 , Aho for 


some datasets (G05; S07), while we plotted only the 
A 5100 > ergs s“^ data for the rest of references that 

meet <10% in host contamination (Sll). Meanwhile, 
the broad Lho could contain the narrow component for 
AKARI data, while the broad and narrow line luminosi¬ 
ties are combined for the G05, S07, and Sll data too. We 
find that the contribution from the narrow component 
to Lfja estimated from section 3.1 and the references, is 
negligible (2% and < 10%, respectively), allowing us to 
consider Ajja to be approximately the line luminosity of 
the broad line component. 


Figure 11 shows a remarkable correlation between 
L 5100 and Aho even when AGNs are drawn from various 
samples covering a wide range of redshifts and luminosi¬ 
ties. This strongly suggests that the physics governing 
the correlation is the same for low and high luminos¬ 
ity AGNs, and there is no strong evolution in the re¬ 
lation from z = 0 to z = 6 over the range explored here. 
The AKARI data points are mildly below the relation 
where the offset could indicate a growing po pulation of 
weak emission line quasars at hig h redshift (jFan et al.l 
119991: iDiamond-Stanic et al.ll200^ . but the overall de¬ 
viations are within the scatter of the data. The possi¬ 
ble downturn in the correlation at very high luminos¬ 
ity (A 5100 lO^'^ergss”^) produces a ^ 0 . 2 dex offset in 
log (Ana/Asioo), but this downturn affects the Mbh es¬ 
timates only by ~0.1dex. The possible downturn could 
be caused by cold accre tion disks of slowly spinning, ex¬ 
trem ely high mass BHs (|Laor fc Davisl[201lt iWang et al.l 
l2014f) . but we will leave the investigation of the possi¬ 
ble downturn as a subject of a future work as its effect 
on Mbh estimates is small. We fitted the relation using 
the linear regression wit h bivariate correlated error s and 
intrinsic scatter (BGES. lAkritas fc Bershad(lll99(TfF^ . to 
find the following result 0 

log AHa ,42 = (1.044 ± 0.008) log A 5100.44 

-f(0.646± 0.011). ^ ^ 

The best-fit relation fits the entire data with 
Cint = 0.095 dex. The flux-flux relation of 5100 A con¬ 
tinuum and Ha shows almost identical slope and the in¬ 
trinsic scatter, suggesting that the tight correlation in 
Equation (2) is not due to a sample selection effect. The 
data points at three redshift intervals, namely 0 < z < 
0.8 (G05, S07, and Sll) for 41.7 < logAgioo < 46.4, 

Throughout this paper we use the BCES fit to derive the 
linear relations. 

Throughout this paper we use subscript numbers to 
the luminosity to indicate its wavelength and unit, such as 
I'SIOO .44 = 1-5100 a/IO'*'* ergss“h 
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Fig. 13.— The FWHMn^-FWHMHa) FWHMH^-FWHMMgn ? and FWHMh^-FWHMciv relations of AGNs from combined references. 
The references abbreviated on the plot are summarized in Table 4. The AKARI data are not present since the sample do not have 
simultaneous coverage of the FWHMs plotted. We limited the G05 data to L 5100 < 10"^^ ®^ ergss”^ to avoid overlap with the Sll data. We 
assign a 10% error for FWHMjj^, FWHMcj^.^ of the N07 data, following their argument. We removed one object in S07 data that overlaps 
with WOO, and two objects in S07 that overlap with P13. The FWHMh/?—FWHMhck relation of G05, and FWHMjj^-FWHMMgji relation 
of WOO are shown in dotted lines, while the relations from this work are shown as solid lines. The zeropoint offset and rms scatter of each 
literature data with respect to our relations, are indicated. The FWHMs are of the broad emission line. 


0.7 < z < 2.2 for 44.8 < logLsioo < 47.0 (H12, S12, 
and M13), and 3.3 < z < 6.2 for 46.0 < logLsioo < 47.2 
(J15), overlap with each other and show no evolution. 
This suggests that Equation (2) is universal, and not 
due to a distance effect like the Malmquist bias. 

To examine the universality of the relation further, we 
discuss how the L^ioQ-Lua relation of G05 at z ~ 0 of 
lower luminosity AGNs and S12 at z^2 of higher lu¬ 
minosity AGNs fare with each literature values. First, 
the G05 relation can describe the Lsioo-Tho relation 
of S07, Sll, and M13 AGNs over the overlapping lu¬ 
minosity interval (Tsioo ^ 10^® ergss“^). When extrapo¬ 
lated to higher luminosity, it starts to deviate from the 
data points regardless of redshift. Likewise, the S12 
relation can describe the Lsioo^Aho, relation down to 
.^5100 ergs s“^ including the S07 data at z<0.4 

and our AGNs at z > 3.3. However, when applied to 
the entire datasets, the G05 and S12 relations show de¬ 
viations from the data at high and low luminosity re¬ 
gions, respectively. Consequently, both relations produce 
Cint = 0.15-0.24 dex against the data which is worse than 
0.095 dex of our L 5 ioo“.^Ha relation. 

To check if the inconsistency in the G05 and S12 
relations at the faint and luminous end arises from a 
possible break in the relation itself, we considered the 
case where the slope changes at L 5100 ~ 10 ^"^'® ergss“^, 
where the G05 and S12 relations meet. For this, 
we computed the CTint of the L 5100 < lO"'^'^ ergss“^ and 
-^5100 > 10"‘'*'^ ergss“^ data against the G05, S12, and our 
relation in Equation (2), and examined if our simple re¬ 
lation is any worse than the combination of G05 and S12 
relations with a break at L 5100 = 10'*'^'^ ergss“^. As in¬ 
dicated in the lower panel of Figure 11, the Uint values 
against our relation is comparable to or slightly smaller 
than the Uint against the G05 or S12 relations at low and 
high luminosities respectively. This suggests that there 
is no strong need for a broken power-law form of the 
Asioo^ina relation, and a simple relation of Equation 
( 2 ) can be employed to describe the response of broad 
line region to the incident continuum emission over the 
covered redshift and luminosity ranges in Figure 11 and 
Table 4. 

Likewise, we plot in Figure 12 the Lsioo^^isso and 
A 51 oo “^3000 relations from the luminosities derived in 
section 3 and taken from references. Apart from the lit¬ 


erature data where the host contamination in L 5100 is 
estimated to be <20% (S07), or minimized by Hubble 
Space Telescope (HST) observations, we limit the litera¬ 
ture sample with L 5100 > ergss“^ to keep the host 

galaxy contaminatio n below 10% ISll) . Me anwhile, we 
replac ed the L 5100 of lWang et al.l (1200911 andiPark et ^ 
(j2r)13[l with the HST data from iBentz et al.l (j2013f l. while 
only including the AGNs with <20% host contamination 
in L 5100 . Similarly to the Lsioo^-^^Hq relation, we do not 
find any evolution in the ^ 5100-^1350 and ^ 5100-^3000 
relations for a particular set of data, and we find the 
best-fit correlation to be, 

logLi350,44 = (0.974 ± 0.023) logL5ioo,44 
±(0.391 ±0.053) 
log ± 3000,44 = (0.973 ± 0.010) log ± 5100,44 
±(0.287 ±0.013). 

The CTint values (dex) to these best-fit relation, as well as 
the relation with respect to S12 are presented in Figure 
12. Like for the ±Ha-± 5 ioo relation, the flux-flux re¬ 
lations show virtually identical slopes and intrinsic scat¬ 
ter to the luminosity-luminosity relations, showing again 
that sample selection effect is not a main driver for the 
relations in Equation (3). 

Finally, we compare the broad line FWHM of H/3, Ha, 
Mg II, and G IV in Figure 13 in order to calibrate the 
Mbh from multiple line based recipes and to check for 
any evolution in the FWHM relations. For good com¬ 
parison of FWHMs, we restricted the mixed samples to 
have the fractional errors of FWHM less than 20%, while 
additionally limiting the S/N and reduced chi-square of 
the SDSS spectra to be S/N>20 and xl < 2- In Figure 
13a, we fit the FWHMn/j-FWHMHa relation from the 
collected data, where we find the offset of each reference 
data to this relation to fall within each scatter. The CTint 
of all the data to our relation, 0.063 dex, is smaller than 
when using the relation from G05. 

Second, we derive the FWHMn^g-FWHMMgn relation. 
Since it is debatable whether to subtract the narrow com¬ 
ponent for t he Mg II line wid th measurement (e.g., Sll), 
we followed lJun fc Iml (|2013f) to average the FWHMMg„ 
derived with and without the subtraction of the narrow 
component. Also, the slope and constant of the relation 
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Fig. 14.— The comparison of Balmer Mbh relations of AGNs, from combined references. The references abbreviated on the plot 
are summarized in Table 4. All measurements are converted to our adopted cosmology and /-factor. We limited the G05 data to 
^5100 < ergss“^ and the Sll data to L 5100 > ergss“^, in order to avoid overlap. The linear fit to the Mbh relation and a 

1-1 relation are represented by solid and dotted lines respectively, while the zeropoint offset between the masses (Azpt), intrinsic scatter 
with respect to a 1-1 relation (cri^t, dex), and the fraction where the masses overlap within error (pcon) are printed. 


are consistent within uncertainty with those derived with 
and without the subtraction of the narrow FWHMMgn 
component. Therefore, we combined the literature data 
with respective treatment of the narrow Mg II and the 
averaged FWHMMgn in SII, to derive the FWHMh,s- 
FWHMMgn relation altogether. We note that the rela¬ 
tive FWHM offset of the literature data to our relation 
shown in Figure 13b are within the scatter of data points, 
indicating that the det ails of fitting to exclude the narrow 
Mg I I component (S07: lMcGill et al.ll2008HDietrich et alJ 
1200911 or to include but subtract it ( Wang et al.l 120011 
S12) do not affect the line widths significantly when com¬ 
pared overall. The rms of all t he data to our rel ation, 
0.085 dex, is similar to that from lWang et al.l (j2009ll . and 
small enough to regard the FWHMMgn as a marginally 
good substitute of FWHMh /3 as much as FWHMhq ■ 

Third, we derive the FWHMh/3-FWHMciv relation in 
Figure 13c. The data can be fitted altogether with a log- 
linear relation, but the CTint of the data to the relation, 
0.212 dex, is large and comparable to the systematic un¬ 
certainty of single-epoch Mbh estimators when scaled as 
the FWHM squared. We checked the effect of fitting 
methodology by deriving all the L-L, FWHM-FWHM 
relati ons with the FITEXY method (| Tremaine et al.l 
I2002H . to find that the slope of the FWHMh,s-FWHMciv 
relation, 1.798 ± 0.026, shows a meaningfully large differ¬ 
ence to the BCES results. Still, we keep the BCES slope 
since it brings the Balmer and C IV Mbh estimates more 
consistent (section 4.3), and we note that the large CTint 
between the EWHMh,s-FWHMciv relation represents a 
poor correlation at best. 

Overall, we find the relations of FWHMhq, 
FWHMMgn and FWHMciv against the FWHMh/ 3 
without any noticeable evolution for the samples con¬ 
sidered, covering a wide range of luminosity or redshift. 
Therefore, although the data for calibration is missing 
at we use our FWHM relations to calibrate the 

Mbh estimators for general usage in terms of redshift. 


Our derived FWHM relations are as follows FI 
log FWHMh/3,3 

= (1.061 ± 0.013) logFWHMHa.3 + (0.055 ± 0.008) 

= (1.226 ± 0.032) logFWHMMgn.3 - (0.078 ± 0.02iy^ 
= (1.054 ± 0.057) logFWHMciv.3 - (0.024 ± 0.045). 

We note the possibility that the FWHM relations estab¬ 
lished with spectra taken in different epochs could be 
affected by variability of the emission line shape, espe¬ 
cially when the dynamic range of the probed FWHMs are 
narrower than that of the continuum luminosities. Still, 
the offset or scatter of the references with and without 
simultaneously acquired data in Figure 13 do not dif¬ 
fer with each other. Also, we find the average and rms 
scatter in the ratio of the SDSS-HI BOSS over SDSS- 
I/SDSS-H FWHMciv from 12, S/N > 15 objects matched 
with the AKARI sample to be 0.99 ± 0.08, which implies 
the effect of variability to the shift or broadening of the 
FWHMh/ 3 -FWHMciv relation is negligible compared to 
the CTint of the relation by an order of magnitude. 

4.3. Updated Mbh estimators 

Using the scaling relations obtained above, we now 
present an updated set of Mbh estimators based on var¬ 
ious lines. Since the reverberation mapping of H/3 and 
U 5100 forms the basis of mass estimation for AGNs, we 
start from the Mbh estimator that uses the 5100 A lumi¬ 
nosity and H/3 line width to take the following form and 
derive other estimators using scaling relations discussed 
above: 

Mbh = ^i7BLR(A5ioo)( ™^^^ )'. (5) 

Previous secondary calibrations to the Mbh were per¬ 
formed through replacing the (T 5100 , FWHMh/?) of the 
local reverberation mappe d sample or SPSS A GNs, by 
(Lna, FWHMhc) (e.g., IGreene fc Hoi [200H) . (L 3000 , 
FWHMMgii) fe.g.. IMcLure fc JarvisI 120021), or 
FWHMciv) (e.g., iVestergaard fc Peterson! 1200611 . Us¬ 
ing the latest R-L relation from iBentz et al.l (l2013rF^ 


Throughout this paper we use subscript numbers to 
the line width to indicate its unit, such as FWHMh/ 3,3 = 

We shifted the relation to Hubble parameter of Hq = 
70 kms“^ Mpc”^ which is our adopted value throughout this work. 
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Fig. 15.— The comparison of Balmer Mbh t-o the UV Mbh of AGNs, from c ombin ed r eferences. The references abbrev iated on the plot are 
summarized in Table 4, with M04, V06 estimators from lMcLure Dunlod {2003) and I Vestergaard &: PetersonI II2006I') . All measurements 
are converted to our adopted cosmology and /—factor. We use the Mbh(^5100! FWHMhck) for th® Balmer masses, unless when the Ha line 
was uncovered and the Mbh(^5100? FWHMh^) was used. The rest of the figure format follows that of Figure 14. We removed the objects 
unused in Figures 13 and 14. 


and the consta nt for the mass e quation (/-factor, / = 
5.1 ± 1.3) from lWoo et al.l (l2013f l. we derive first the H/3 
Mbh estimator as, 

i?BLR = (34.7 ± 2.5) L 5100.44 (0-533±0.034) 

MbhC^sioo, FWHMh/ 3) = (8.63 ± 2.29) x 10® (6) 

X L5100,44 FWHMh/3.3 "Mq. 

Replacing the H/3 line width through the FWHMh/?- 
FWHMhq, relation from Equation (4) we get the Mbh 
from (L 5100 , FWHMna), 

MBH(i5ioo,FWHMHa) = (1.11 ± 0.30) X 10^ 

X L5100.44 FWHMHa.3 

Finally, the application of Equation (2) yields the Mbh 
from (Lho, FWHMna), 

MBH(FHa, FWHMna) = (5.20 ± 1.41) x 10® 

X LHa.42 FWHMHa.3 

Likewise, replacing the continuum luminosity and line 
width in Equation (6) to those of Mg 11 and C IV using 
Equations (3) and (4) yields, 

MBH(L30oo,FWHMMg„) = (4.19 ± 1.19) x 10® 

X 7.3000.44 ^° ®^®^° °®®^ FWHMMg„.3^"-^''^°-°®^M© 

Mbh(Li350, FWHMcv) = (4.72 ± 1.63) x 10® 

X Li350,44 (0 ®47±0.037) FWHMc.v.3 

The Mbh’s based on Ha and C IV from our estima¬ 
tors for the AKARI sample are given in Table 3. It 


is worth noting that Equation (5) may not hold if the 
FWHMh /3 is not exactly propor t ional to the velocit y dis¬ 
persion, a (|Peterson et ai.l[200l : ICollin et al.ll2006ll . or if 
the R-L relation breaks down at high luminosity where 
the relation has not been tested extensively with the re¬ 
verbe ration technique (L 5100 > 10"^®ergss~^. lBentz et all 
nnn). Nevertheless, considering the advantages of using 
the FWHM (to be robust under poor sensitivity, wings 
in the line profile, or deblending, see S12) for the single 
epoch mass estimation, and the current expectations in 
the high luminosity R-L relation (see section 5.1), our set 
of calibrations have the merit where the rest-UV to op¬ 
tical Mbh estimations are mutually consistent through a 
wide range of redshifts, luminosities, and fitting method¬ 
ologies. 

To check the consistency between our Mbh estima¬ 
tors, we first compare the Balmer Mbh using the es¬ 
timator in this work and existing estimators in Figure 
14. The Mbh,H /3 and Mbh,H a are compared, where a 
AMbh < 0.3 dex error cut is applied to the data out 
of various references. The Mbh values are derived from 
the G05, Sll estimators and from this work (Equations 
(6), (7)) in Figures 14a-14c respectively, where the ze- 
ropoint offset between the Ha and AjS masses (A^pt), 
intrinsic scatter with respect to a 1-1 relation ((Tint), and 
the fraction of the data points where 1 a errors of Mbh 
values from different estimators overlap with each other 
(Pcon) are shown. We find that the masses from our es¬ 
timators are closer to a 1-1 relation than of G05 or Sll, 
throughout the range 10®“^®Mq. Also, fully considering 
the propagated errors in the mass equation, our estima¬ 
tors may overestimate the Mbh uncertainty when com¬ 
paring the H/3 and Ha masses as reflected from the nega¬ 
tive (Tint, which suggests that the actual difference in the 
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TABLE 5 

isiOO-i'Ha RELATION 


Reference 

(1) 

OL 

(2) 

d 

(3) 

Cint 

(4) 

(5) 

N 

(6) 

z 

(7) 

log 7-5100 
(8) 

Dataset 

(9) 

Greene &; Ho f2005) 

0.720 ±0.002 

1.157±0.005 

0.078 

0.241 

229 

<0.35 

41.7-45.0 

G05 

Shen & Liu f2012i 

0.791 ±0.093 

1.010 ±0.042 

0.088 

0.153 

60 

1.5-2.2 

45.4-47.0 

S12 

This work 

0.646 ±0.011 

1.044 ±0.008 

0.095 

0.095 

464 

0.0-6.2 

41.7-47.2 

GU5, yU7, yil, LL12, 
S12, M13, J15 


Note. — The L 5100 —Thq relations from previous studies and this work, where a and /3 are defined as logLHQ ;,42 = o + ^ log L 5100 , 44 - 
Column 1: Reference; Column 2: a and its uncertainty (Itr); Column 3: /3 and its uncertainty; Column 4: Intrinsic scatter (dex) of the 
relation over the L 5100 range covered; Column 5: Intrinsic scatter of the relation over the range 41.7 < logLsioo < 47.2; Column 6: Number 
of objects used; Column 7: Redshift range; Column 8: Range of L 5100 (e^gss”^); Column 9: References of AGN dataset with abbreviations 
from Table 4. 


TABLE 6 

ViRiAL Mass Estimators 


Reference 

(1) 

OL 

(2) 

d 

(3) 

7 

(4) 

N 

(5) 

z 

(6) 

log Mbh 

(7) 

Method 

(8) 

A7bh(A5100,44, FWHMh/3,3) 

Greene & Ho (2005) 

6.64 ±0.02 

0.64 ±0.02 

2 

35 

<0.37 

5.5-9.0 

RM 

Vestereaard Peterson (2006) 

6.91 ±0.02 

0.5 

2 

25 

0.00-0.29 

7.2-9.3 

RM 

Bentz et al. (2013). this work 

6.94 ±0.12 

0.533 ±0.034 

2 

41 

0.00-0.29 

6.0-10.7 

RM 

A7bh(7-5ioo,44, FWHMhq.s) 

Greene Ho [2005) 

6.70 ±0.06 

0.64 ±0.02 

2.06 ±0.06 

162 

<0.37 

5.2-9.0 

SE 

J^hen kr. Liu (*2012') 

7.01 

0.555 

1.87 

60 

1.5-2.2 

8.8-10.3 

SE 

This work 

7.05±0.12 

0.533 ±0.034 

2.12 ±0.03 

654 

0.0-2.4 

5.8-10.6 

SE 

MBH(7Hc.,42,FWHMHa,3) 

Greene Ho [2005) 

6.30 ±0.08 

0.55 ±0.02 

2.06 ±0.06 

243 

<0.37 

5.1-9.0 

SE 

Shen Liu (*2012) 

6.55 

0.564 

1.82 

60 

1.5-2.2 

8.8-10.2 

SE 

This work 

6.72±0.12 

0.511 ±0.033 

2.12 ±0.03 

969 

0.0-6.2 

5.6-10.5 

SE 

A7bh (±3000,44, FWHMMgn, 3 ) 

McLure & Dunloo 12004) 

6.5 

0.62 

2 

22 

0.03-0.37 

7.5-8.9 

RM 

Vesterffaard &; Osmer (2009) 

6.86 

0.5 

2 

- 

- 

- 

SE 

Wanpj et al. (2009') 

7.15 ±0.27 

0.46 ±0.08 

1.48 ±0.49 

29 

0.00-0.16 

6.3-9.0 

RM 


6.95 

0.584 

1.71 

60 

1.5-2.2 

8.8-10.1 

SE 

This work 

6.62±0.12 

0.548 ±0.035 

2.45 ±0.06 

1010 

0.0-6.4 

6.8-10.3 

SE 


M'bh(^1350,44, EWHMciv,3) 


Vestergaard Peterson (2006) 

6.66 ±0.01 

0.53 

2 

27 

0.00-0.23 

5.5-9.3 

RM 

Shen & Liu [2012) 

8.02 

0.471 

0.24 

60 

1.5-2.2 

8.9-9.6 

SE 

Park et 

7.48 ±0.24 

0.52 ±0.09 

0.56 ±0.48 

25 

0.01-0.23 

7.0-9.0 

RM 

This work 

6.67±0.15 

0.547 ±0.037 

2.11±0.11 

258 

0.0-5.4 

6.8-10.7 

SE 


Note. — The Mbh relations from previous studies and this work, where o, 7 are defined as log Mbh = a + log L + 7 logFWHM. 
Column 1: Reference; Column 2: oc and its uncertainty (Icr); Column 3: ^ and its uncertainty; Column 4: Column 3: 7 and its uncertainty; 
Column 5: Number of objects used for luminosity or line width calibration; Column 6: Redshift range of objects used for luminosity or line 
width calibration; Column 7: Range of Mbh from reference; Column 8: Method of calibration, where RM, SE denote calibrations of the 
luminosity and line width using the reverberation mapped and single epoch samples, respectively. 


/-factor or -Rblr between the H/3 and Ha line emitting 
regions are likely to be smaller than their uncertainty. 
Therefore, we regard the H/? and Ha Mbh from our esti¬ 
mators to be indistinguishable, which is supported by the 
high fraction (98%) of H/3/Ha masses to be consistent 
within measurement error. 

With the H/3 and Ha Mbh recipes checked to be mu¬ 
tually consistent, we further compare the Balmer to the 
UV-based Mbh with AMbh < 0.3 dex, in Figure 15. 
We use Mbh.H a as the Balmer mass due to the stronger 
Ha emission than H/3 and the availability of Ha line in 
our AKARI data, unless only the Mbh,h /3 can be esti¬ 
mated using the measurements in the literature. Com¬ 
pared to the conventional estimator s calibrated at rel¬ 
ative l y low luminosity and redshif t (jMcLure fc Dunlcnil 
12004 IVestergaard &: Petersonll200^ or relatively high lu¬ 


minosity and redshift (S12), our calibrations (Equations 
(6), (7), (9), (10)) bring the rest-UV (Mg 11 in Figures 
15a-i5c, C IV in Figures 15d-15f) and Balmer Mbh’s 
to be closer to a 1-1 relation in the range in 

terms of both zeropoint offset and intrinsic scatter. Also, 
we find that our Mbh estimators fully considering the 
propagated errors give more realistic value of the error 
bars compared to estimators that do not, and raise the 
Peon- Moreover, previous UV-based Mbh’s are tilted to 
ours largely from the different slope of the FWHM term 
to that of our estimato r (e.g., 0.24 in S 12. 2.11 in this 
work for C IV, see also, iPark et Tl 120131 1. and this pro¬ 
duces a systematic offset in Mbh values. For example, 
we find that S12 estimators produce a systematic differ¬ 
ence between C IV and Balmer Mbh in such a way that 
Mbh,Civ are larger by 0.7 dex at MBH,Ha/H /3 = lO'^M© 
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and smaller by 0.4 dex at MBH,Ha/H /3 = (Figure 

15e). Previous UV-based Mbh calibrations from lim¬ 
ited dynamic range yields a-mt to the Balmer Mbh as 
small as our estimator, and may place the C IV and 
Balmer masses to be consistent within a narrow range 
of masses, but we caution on the usage of conventional 
estimators when extensively comparing the rest-optical 
and UV Mbh- 

Even when using our rest-UV mass equations with 
minimized zeropoint and systematic offsets from a large 
dynamic range, the intrinsic scatter is another issue. Al¬ 
though our (Tint of the MBH,Mgn to the Balmer Mbh is 
relatively small (0.09dex), the (Tint of the Mbh.Civ to 
the Balmer Mbh (0.40 dex) is comparable to the system- 
atic uncertainty of the single epoch mass estimator itself 
(|Bentz et al.ll2013l : lWoo et al.ll2013ll . Moreover, since our 
estimators take into account the measurement error of /- 
factor, R — L relation, scatter in the UV-optical luminos¬ 
ity and line width relations, the non-negligible fraction 
of rest-UV masses that deviate from the Balmer masses 
more than their overlapping measurement errors (13% of 
the MBH.Mgn and 41% of the Mbh.Civ) suggest that one 
needs to be cautious about the usage of the C IV line 
when deriving Mbh values. 

5. DISCUSSION 

5.1. Reliability of single-epoch AGN black hole masses 

Does the relations between continuum and line lu¬ 
minosities, and the FWHM relations that hold up 
to the highest luminosity AGNs, guarantee the ac¬ 
curate mass estimation of the most luminous AGNs? 
Unfortunately the answer is not clear yet, since the 
-Rblr.h/sMsioo relation is not observationally probed for 
-^5100 > 10"^® ergs s“^ (jBentz et al.ll2013ll . A hint to es¬ 
timate the high luminosity end behavior of the opti¬ 
cal R-L relation is to look at the i?BLR,Civ“-^i 350 re¬ 
lation, where the i?BLR,Civ traces the inner part to the 
H/3 line region. Although the curren t number of G IV 
reverberation measurements is small (iKasoi et al]l2007L 
iSluse et all 1201 IL iGhelouche et al.l 1201211 . th^ cover up 
to luminosities of L 1350 = 1047.O 

ergss ^ and £.-sinn = 
1046.8 ergss~^. The slope of the relation in iKaspi et al.l 
(|2007ll . (0.52-0.55)±0.04, is within the un certainty to 
the s lope of the i?BLR,H/ 3 “A 5 ioo relation (|Bentz et al.l 
1201311 . 0.533 ± 0.034. This suggests that the ratio of 
-RBLR,H/3-to-i?BLR,Civ is nearly a constant at 2-4 over 
A 5100 < 10"*® ergss“^, where the proportionality of T 5100 
and Ti 35 o suggests that it is likely so at a higher lumi¬ 
nosity range. Therefore, we do not expect a break in 
the luminous end optical R-L relation unless the optical 
data will somehow fail to form a simple relation at high 
luminosity where the current UV relation holds. Future 
compilations of both the G IV and H/3 broad line region 
sizes will help to better constrain the luminous end R-L 
relation. 

Even if we find it plausible to assume that the 
??blr,h,S“T 51 oo relation does not change its slope in the 
luminous end, the large scatter in between the Balmer 
and rest-UV Mbh imposes further limitations on the us¬ 
age of UV Mbh estimators. As discussed in the introduc¬ 
tion, it has been known that there is a large (Tint in the 
G IV and Balmer Mbh relation in comparison to that 
of the Mg II and Balmer. Our comparison of Mbh, Ha 


against Mbh,Civ high luminosity AGNs at z > 3.3 
further shows consistent results. Indeed, we find 4 out of 
11 Mbh, Civ of AKARI quasars in Figure 15f to be scat¬ 
tered to the Mbh, Ha by an order of magnitude or larger, 
therefore the use of G IV masses can be questioned up 
to the highest redshifts even when considering the ad¬ 
ditional uncertainties in the AKARI Mbh, Ha tested in 
Figure 5. Overall, our result supports that the Mbh es¬ 
timate from G IV contains a rather large uncertainty at 
any redshift or luminosity. 

The greatest uncertainties when calibrating the UV 
Mbh estimators in this work ((Tint > 0.1 dex), come from 
the scatter in the T5ioo“-^i350 and FWHMh/ 3 -FWHMciv 
relations. Since a variety of obscuration in the rest-UV 
continuum of quasars may result in a scattered ^ 5100 “ 
Li 350 relation, we suggest to take into account the level 
of obscuration when establishing the Mbh,Civ estimator, 
although careful treatment is required as it is controver¬ 
sial whether the color correction into the Mbh works to 
reduce the scatter between the G IV and Balmer masses 
l|Assef et al.l 120111 : S12). Meanwhile, the scatter in the 
FWHM relations which is the largest for the FWHMh^- 
FWHMciv relation, could arise from broad absorption 
features or a no n-reverberating component of the C IV 
broad emission (iDennevI 1201211 . where future reverber¬ 
ation mapping will further enable us to test. Without 
detailed analysis on the origin of the scatter in between 
the Balmer and rest-UV Mbh, we caution on the indi¬ 
vidual measurement of Mbh from UV, especially the one 
using G IV. 

5.2. On the massive end black hole mass evolution 

Having cross-calibrated the Mbh estimators, we plot 
the Balmer and Mg II Mbh measurements of quasars 
along redshift in Figure 16a, from our AKARI obser¬ 
vations and the literature indicated on the figure. First 
of all, we notice that the most massive envelope that 
stays at 10 ^°Mq up to z ^ 5, starts to disappear at 
higher redshift. To quantify how significant the Mbh’s 
are evolving at the massive end for z > 5 quasars, we 
performed the Kolmogorov-Smirnov (K-S) test to com¬ 
pare the mass distribution at above and below z = 5, at 
Mbh > 2 X 10^Mq. This mass threshold for complete¬ 
ness is governed by the shallow AKARI data, through 
the L 5100 ^ 10 "^®'® ergss“^ limit at S/NHa = 3 and the 
FWHM limit of 2500kms“^ from section 3.1. The K- 
S probability pxs = 0.34% turns out to be meaning¬ 
fully small, supporting the massive end BHs at z < 5 to 
be heavier than at z > 5. We checked our results to be 
unchanged much by the inclusion of AKARI data that 
could bias the masses by up to ^0.1 dex (section 3.1), 
by performing the same test while giving a — 0.1 dex cor¬ 
rection to the AKARI data, only to find the trend to be 
mildly weaker {pxs = 0.64%). 

We cross-check if this evolution in Mbh is also re¬ 
flected through the Mbh histograms. In order to cor¬ 
rect for the luminosity selection effect in number count¬ 
ing, we plot the 1450 A absolute magnitude (M1450) his¬ 
togram in Figure 16b. We overplot the quasar lumi¬ 
nosity functions in redshift bins of2<z<3, 3<z<5, 
z > 5, scaled by a constant so that they match the M 1450 
counts in total number at M 1450 < —26.4. We use only 
the Mi 45 o < —26.4 data for better completeness in seal- 
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Fig. 16.— The Mbh of AGNs along redshift (left). The Mbh from Ha are in circles, H/3 in squares, and Mg II in triangles, while our 
AKARI data points are in red circles. For objects with multiple lines, we choose the Ha over H/3, and H/3 over Mg II, while for repeated 
measurements of an object we used the latest results. Most of the refe r ences abbreviated o n the plot are summarized i n Table 4, with 
S04, K07, WIO, Dll, Til, D14, W15 additionally from IShernmer_et_al.l j200J) , IKurk et al.l II2007I1 , IWillott et al.l Il2010al l , IDe Rosa et al.l 
II2011I1 . [Trakhtenbrot et al.l 11201111 . FPe Rosa et al.l ||2014I1 . and IWu eT^TTKoT^ . We do not plot the error for the S04 data as they are not 
available, and do not use a part of the Til data with poor quality (flagged “3”). The model tracks of exponential BH growth that matches 
the massive envelope of AGNs at 2 ~ 5 and 7 are plotted, with each set of parameters describing the curve shown. On the right, we plot 
the histogram of M 1450 (top) and Mbh (bottom) at each denoted redshift bin in circles, with errors determined through Poisson statistics. 
The luminosity functions scaled to the observed number counts to yield the same number at M 1450 < —26.4 at each redshift (top), and the 
resultant mass counts from the luminosity scaling and a FWHM > 2500 kms“^ cut (bottom), are overplotted in lines with styles marked 
next to the redshift. 


ing the luminosities. The break luminosity (M 1450 ), the 
faint end slope (a), and the bright end slope (/3) of the 
luminosity function are fixed to M 1450 = —26.39 at z > 3 
and Mi 45 o = —2 5.5 at z < 3 while (a, / ?)=(—1.80, —3.26 ) 
at all redshifts (|McGreer et alJ [20131 : iRoss et ^120131) . 
The set of (a, /3) are fixed as the evolution in the a value 
does not affect the luminosity regime of interest much, 
and P stays nearly constant with redshift. Overall, we 
find that the quasar samples at z < 5 are in excess of 
luminous objects compared to the luminosity functions. 

Not only there is a luminosity selection effect in the 
quasars plotted in Figure 16a, there is a FWHM limit for 
each reference data which is the poorest for the AKARI 
data. Taking into account of these effects, we plot the 
Mbh histograms in Figure 16c before (points) and af¬ 
ter (histograms) applying the luminosity selection cor¬ 
rection and a FWHM = 2500 km s“^ cut. After correct¬ 
ing the luminosity selection and giving the FWHM cut 
the Mbh distribution is shifted to more IO^Mq BHs and 
less IO^^Mq BHs at z < 5. Still, the number of IO^^Mq 
BHs start to drop at z > 5, whereas the IO^Mq BHs do 
not exhibit an evolution up to Poisson error. Thus, we 
interpret the Mbh histogram that there is an evolution in 
the number of 10^°Mq BHs to be increasing with time 
at z > 5. The z > 5, ~ IO^Mq AGNs are likely to be 
in a rapidly growing state such that their masses can 
soon reach ~ 1 O^°M 0 , consistent with the higher Ed- 
dington ratio trend f ound in distant, lum inous AGNs 
(|Willott et al.l [ 2010 ^ IDe Rosa et al.l[ 2 ?)TTI ). Note that 
the number density of quasars evolve strongly from z = 2 


to 6 (e.g., IWillott et al.ll20l'0^ . We caution that the 
mass and luminosity histograms in Figures 16b and 16c 
do not reflect this number density evolution, and that 
Figures 16b and 16c are used for comparing the shape 
of these functions. The fact that the numbers of z = 2 
through 6 quasars appear similar in these figures is due 
to the fact that the sampling of quasars is different at 
each redshift bin. 

We further investigate the upper mass envelope of z > 5 
AGNs to check if their extremely high masses can be 
explained by the earliest BH growth from mode l seed 
masses. For this we follow Iwionteri fc R.eesI (|2005D to as¬ 
sume a continuous, Eddington limited accretion of mat¬ 
ter to a BH seed at z = 20, where the BH will grow with 
time as Equation (1), with to = tz= 2 o and Mq = Mz= 2 o 
for the age of the universe and the seed BH mass, re¬ 
spectively at z = 20. Eor each given e = 0.05,0.1,0.2, 
we determined a pair of seed masses Mq such that the 
growth curve reaches the observed massive limit at z ~ 5 
and z ~ 7 respectively, shown in Figure 16a. 

When e = 0.05, the seed masses to explain the most 
massive BHs at z ^ 5 and 7 are Mg ~ and 

1 O“^M 0 respectively, which in other words, gives plenty 
of time for seed BHs of any meaningful mass to grow up 
to the most massive quasars without a strictly continu¬ 
ous, Eddington limited accretion. One can expect to find 
fully grown massive BHs at z ~ 6 detectable as high or 
low Eddington ratio AGNs, if such rapid growth is pos¬ 
sible. Luminous AGNs at z ^ 6 have high Eddington ra¬ 
tios so far, but future studies may uncover fainter objects 
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with weakly active BHs. If e = 0.1 that is rough ly con- 
sistent with the Soltan argument measurements (|Soltanl 
119821 : 1^1 fc Tremainell200^ . it predicts Mq 
which gives a reasonable estimate of the seed mass to 
be Population III stars. However, it may be difficult 
for stellar seed BHs to keep its Eddington limited accre¬ 
tion from z = 20 to 7-5, which is 0.6-1.0 Gyrs in duration 
and is longer than a typical quasar lifetime (<0.1Gyr, 
iHopkins et al.l[^005f) . perhaps requiring extended period 
of accretion or BH-BH mergers. Lastly, the e = 0.2 
model which corres ponds to the case when BHs a re ro¬ 
tating rapidly le.g.. lKeri1ll963tlMarconi et al.ll2004[l only 
accepts very heavy seed BHs at Mq ~ 1O®-’^“®'®M0, re- 
quiring supercritical acc retion from lighter seed masses 
(iVolonteri fc Reed 1200511 or supports direct co llapse of 
massive primordial gas (|Volonteri et al.l 1200811 . With 
e = 0.2, we expect that the most massive ^lO^*' Mq BHs 
from Mq ~ IO^Mq seeds can only start appearing at the 
redshifts between 5 and 6. 

6 . SUMMARY 

We measured the redshifted Ha emission and rest- 
UV to rest-optical continuum properties of 155 luminous 
quasars at 3 < z < 6 using the AKARI spectra and other 
datasets, and estimated their Mbh- We summarize our 
findings as the following. 

1. The Lsioo^Lhci relation holds up to the most lumi¬ 
nous quasars (L 5100 10"*^ergss“^) with a single slope 

unchanging up to z ~ 6, suggesting a consistent response 
of the broad line region to the incident continuum irre¬ 
spective of AGN luminosity and redshift. 

2. The relations between rest-optical and UV continuum 
and line luminosities, and the FWHM relations hold up 
to the highest luminosity AGNs. Together with predic¬ 
tions of an extended optical R-L relation to the highest 
luminosities, it enables us to calibrate the rest-UV Mbh 
estimators to be consistent with the Balmer masses over¬ 
all. However, some of the rest-UV and optical Mbh are 
scattered more than their uncertainties including the er¬ 
rors from the recipe, for only 13% of Mg II masses but 
for 41% of the G IV estimation. The G IV masses have 
a 0.40 dex intrinsic scatter to the Balmer masses, which 
places negative implications on its reliability. 

3. The massive end envelope of Mbh steeply evolves 

at z > 5, suggesting they are in a rapidly growing state 
from given seed masses. The most massive BHs at z = 5- 
7 can be explained by the Eddington limited, continu¬ 
ous accretion onto ~ seed masses at z = 20 if 

e = 0.1 and the formation redshift is z = 20 for the seed 
BH, while there are a range of viable accretion rates and 
seed masses if e is different. 

We expect future observations to compile sensitive 
measurements of (L 5100 : Lho) to better identify the out¬ 
liers in the relation, where an example would 

be a small population of weak emission line quasars. 
Also, future rest-optical reverberation mapping of high 
luminosity AGNs will verify if the prediction on the R- 
L relation to extend with a single slope will hold. Most 
importantly, the discovery of the highest redshift quasars 
will further uncover the evolutionary tracks of the earli¬ 
est BH growth, also improving the current understanding 
of Mbh growth at z > 5 from small number statistics. 
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project of the University of California, Los Angeles, 
and the Jet Propulsion Laboratory/California Institute 
of Technology, funded by the National Aeronautics and 
Space Administration. 
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TABLE 7 

CONTINUM AND LINE BASED PROPERTIES OF THE SAMPLE 


Name 

(1) 

-^ref 

(2) 

^Ha 

(3) 

log Li350 

(4) 

log ±5100 

(5) 

log Lhc 
(6 ) 

FWHM3,c,v 

(7) 

FWHM3,hc< 

(8) 

log Mbh.Civ 

(9) 

log MbH.Hq: 
(10 ) 

SDSS J000239.39+255034.8 

5.80 

5.79 

99.00 ±99.00 

46.56 

±0.02 

45.14± 

0.11 

99.00 

± 

99.00 

2.50 

±- 

-1.00 

99.00 ±99.00 

9.25 

±- 

-1.00 

Q 0000-26 

4.10 

4.11 

47.40 ±0.01 

47.13 

±0.01 

45.48 ± 

0.13 

4.33 

± 

0.28 

2.50 

±- 

-1.00 

9.88 ±0.22 

9.56 

±- 

-1.00 

SDSS J000552.34-000655.8 

5.85 

5.85 

99.00 ±99.00 

46.04 

±0.08 

44.75 ± 

0.32 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

BR J0006-6208 

4.45 

4.49 

46.90 ±0.02 

46.71 

±0.00 

45.46 ± 

0.06 

11.33 

;± 

:1.33 

2.85 

± 

0.85 

10.48 ±0.24 

9.46 

± 

0.31 

SDSS J001115.23-I-144601.8 

4.97 

99.00 

99.00 ±99.00 

46.78 

±0.02 

47.97 ± 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

BR J0018-3527 

4.15 

4.15 

99.00 ±99.00 

46.72 

±0.01 

45.39 ± 

0.29 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

PMN J0022-0759 

3.90 

99.00 

99.00 ±99.00 

46.57 

±0.01 

47.97 ± 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

BR 0019-1522 

4.53 

99.00 

46.95 ±0.02 

46.58 

±0.02 

47.92 ± 

-1.00 

3.48 

± 

0.33 

99.00 

± 

99.00 

9.43 ±0.21 

99.00 

± 

99.00 

BR J0030-5129 

4.17 

99.00 

99.00 ±99.00 

46.55 

±0.01 

48.02 ± 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

BRI J0048-2442 

4.15 

99.00 

99.00 ±99.00 

99.00: 

±99.00 

48.08 ± 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

BRI J0048-244 

4.15 

99.00 

99.00 ±99.00 

46.07 

±0.03 

47.90 ± 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

SDSS J010619.24-I-004823.3 

4.45 

4.44 

46.99 ±0.01 

46.80 

±0.02 

45.69 ± 

0.19 

3.00 

± 

0.20 

7.75 

± 

2.92 

9.32 ±0.11 

10.42 

;± 

:0.38 

BRI J0113-280 

4.30 

99.00 

99.00 ±99.00 

46.38 

±0.01 

47.92 ±. 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

SDSS J011351.96-093551.1 

3.67 

3.68 

46.93 ±0.03 

46.54 

±0.01 

47.65 ±- 

-1.00 

8.64 

± 

0.36 

99.00 

± 

99.00 

10.25 ±0.22 

99.00 

± 

99.00 

SDSS J011453.81-103934.0 

3.73 

99.00 

46.65 ±0.12 

46.62 

±0.01 

47.69 ± 

-1.00 

5.11 

± 

0.35 

99.00 

± 

99.00 

9.62 ±0.22 

99.00 

± 

99.00 

SDSS J012211.11-1-150914.3 

4.46 

99.00 

47.00 ±0.04 

46.58 

±0.02 

47.91 ± 

-1.00 

7.59 

± 

0.82 

99.00 

± 

99.00 

10.17 ±0.23 

99.00 

± 

99.00 

SDSS J012403.77-I-004432.7 

3.83 

3.83 

47.11 ±0.01 

46.83 

±0.01 

45.44 ± 

0.20 

8.06 

± 

0.29 

3.76 

± 

1.69 

10.29 ±0.15 

9.77 

± 

0.44 

SDSS J012700.68-004559.1 

4.08 

4.07 

47.01 ±0.01 

46.68 

±0.01 

45.14 ± 

0.33 

6.96 

± 

0.27 

99.00 

± 

99.00 

10.10±0.13 

99.00 

± 

99.00 

BRI J0137-422 

3.97 

3.98 

99.00 ±99.00 

46.42 

±0.01 

45.49 ± 

0.22 

99.00 

± 

99.00 

5.95 

± 

2.46 

99.00 ±99.00 

9.98 

± 

0.41 

SDSS J014049.18-083942.5 

3.71 

3.69 

47.24 ±0.01 

46.96 

±0.01 

45.51 ± 

0.12 

4.84 

± 

0.07 

5.05 

± 

1.27 

9.89 ±0.21 

10.11 


:0.28 

SDSS J015048.83-I-004126.2 

3.70 

3.70 

46.95 ±0.01 

46.64 

±0.01 

45.42 ± 

0.17 

4.14 

± 

0.34 

6.47 

± 

2.18 

9.58 ±0.11 

10.17i 

:0.35 

SDSS J021646.94-092107.3 

3.72 

99.00 

47.01 ±0.01 

46.74 

±0.01 

45.47 ± 

0.11 

9.87 

± 

0.96 

99.00 

± 

99.00 

10.42 ±0.18 

99.00 

± 

99.00 

SDSS J023137.65-072854.4 

5.41 

99.00 

46.66 ±0.07 

46.59 

±0.02 

48.20 ± 

-1.00 

3.03 

± 

0.33 

99.00 

± 

99.00 

9.15 ±0.22 

99.00 

± 

99.00 

BR J0234-1806 

4.31 

99.00 

99.00 ±99.00 

46.54 

±0.01 

47.99 ± 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

SDSS J024447.79-081606.0 

4.07 

99.00 

47.19 ±0.02 

46.94 

±0.01 

47.96 ± 

-1.00 

9.76 

± 

0.54 

99.00 

± 

99.00 

10.51 ±0.23 

99.00 

± 

99.00 

BR J0301-5537 

4.13 

99.00 

99.00 ±99.00 

46.41 

±0.01 

47.86 ±- 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

BR J0307-4945 

4.73 

4.78 

47.47 ±0.01 

47.12 

±0.01 

45.79 ± 

0.11 

6.40 

± 

0.86 

7.59 

± 

1.51 

10.27 ±0.25 

10.57 i 

:0.24 

BR J0324-2918 

4.62 

4.60 

99.00 ±99.00 

46.94 

±0.01 

45.29 ± 

0.27 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

SDSS J034402.85-065300.6 

3.96 

99.00 

46.98 ±0.02 

46.44 

±0.01 

47.83 ± 

-1.00 

7.23 

± 

0.24 

99.00 

± 

99.00 

10.11 ±0.21 

99.00 

± 

99.00 

BR J0426-2202 

4.32 

4.33 

47.13 ±0.01 

46.72 

±0.01 

45.46 ± 

0.20 

5.55 

± 

1.16 

4.21 

± 

1.80 

9.96 ±0.28 

9.82 

± 

0.42 

BR J0525-3343 

4.38 

99.00 

47.01 ±0.02 

46.84 

±0.01 

47.89 ±- 

-1.00 

11.23 

;± 

:2.14 

99.00 

± 

99.00 

10.54 ±0.28 

99.00 

± 

99.00 

BR J0529-3526 

4.41 

4.41 

99.00 ±99.00 

46.64 

±0.01 

45.35 ± 

0.33 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

BR J0714-6455 

4.46 

4.44 

47.32 ±0.01 

46.84 

±0.01 

45.45 ± 

0.23 

11.4S 

!± 

:0.33 

7.85 

± 

3.33 

10.72 ±0.23 

10.46 

i± 

:0.42 

SDSS J073149.51-I-285448.7 

3.68 

3.69 

46.96 ±0.02 

46.50 

±0.01 

45.13 ± 

0.28 

7.08 

± 

0.63 

99.00 

± 

99.00 

10.09 ±0.22 

99.00 

± 

99.00 

SDSS J075303.33-I-423130.8 

3.59 

3.56 

47.23 ±0.01 

46.79 

±0.01 

45.32 ± 

0.19 

2.93 

± 

0.08 

6.24 

± 

2.29 

9.42 ±0.10 

10.22 

!± 

:0.37 

SDSS J080430.57-I-542041.1 

3.76 

3.75 

46.99 ±0.01 

46.74 

±0.01 

45.29 ± 

0.24 

7.13 

± 

0.32 

99.00 

± 

99.00 

10.11 ±0.15 

99.00 

± 

99.00 

SDSS J080849.43-I-521515.3 

4.46 

99.00 

47.16 ±0.01 

46.81 

±0.01 

47.91 ± 

-1.00 

8.41 

± 

0.31 

99.00 

± 

99.00 

10.35 ±0.13 

99.00 

± 

99.00 

SDSS J081754.52-I-413225.3 

3.54 

99.00 

46.59 ±0.02 

46.40 

±0.01 

47.67 ± 

-1.00 

4.94 

± 

0.12 

99.00 

± 

99.00 

9.55 ±0.11 

99.00 

± 

99.00 

SDSS J081827.40-I-172251.8 

6.00 

99.00 

99.00 ±99.00 

46.74 

±0.12 

48.40 ±- 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

SDSS J083118.52-I-424728.8 

3.32 

99.00 

46.58 ±0.04 

46.27 

±0.01 

47.58 ±- 

-1.00 

7.21 

± 

0.48 

99.00 

± 

99.00 

9.90 ±0.16 

99.00 

± 

99.00 

SDSS J083700.82-I-350550.2 

3.31 

3.30 

46.87 ±0.02 

46.62 

±0.01 

45.36 ± 

0.19 

4.03 

± 

0.17 

6.42 

± 

2.50 

9.52 ±0.12 

10.16 

i± 

:0.39 

SDSS J083946.22-I-511202.8 

4.39 

4.39 

46.84 ±0.03 

46.71 

±0.01 

45.61 ± 

0.21 

5.78 

± 

0.21 

6.22 

± 

2.50 

9.83 ±0.15 

10.17i 

:0.40 

SDSS J084119.52-I-290504.4 

5.96 

99.00 

99.00 ±99.00 

46.18 

±0.07 

48.35 ± 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

SDSS J085837.95-I-052141.7 

3.53 

99.00 

46.84 ±0.03 

46.53 

±0.01 

47.62 ± 

-1.00 

6.82 

± 

0.68 

99.00 

± 

99.00 

9.98 ±0.22 

99.00 

± 

99.00 

87GB 090153.2-1-694215 

5.47 

99.00 

99.00 ±99.00 

46.26 

±0.03 

48.11 ± 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

SDSS J090634.84-I-023433.8 

4.51 

99.00 

47.03 ±0.02 

46.73 

±0.02 

47.99 ± 

-1.00 

8.95 

± 

1.18 

99.00 

± 

99.00 

10.34 ±0.19 

99.00 

± 

99.00 

SDSS J092721.82-1-200123.7 

5.77 

99.00 

99.00 ±99.00 

46.20 

±0.06 

48.32 ± 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

SDSS J093554.45-I-525616.5 

4.01 

3.99 

46.90 ±0.02 

46.78 

±0.01 

45.43 ± 

0.22 

5.74 

± 

0.37 

3.04 

± 

1.59 

9.86 ±0.21 

9.55 

± 

0.51 

BR 0951-0450 

4.37 

99.00 

46.71 ±0.03 

46.46 

±0.02 

48.03 ± 

-1.00 

5.64 

± 

0.60 

99.00 

± 

99.00 

9.74 ±0.22 

99.00 

± 

99.00 

SDSS J095511.32-1-594030.7 

4.34 

4.30 

47.03 ±0.02 

46.81 

±0.01 

45.72 ± 

0.14 

4.61 

± 

0.26 

3.71 

± 

1.11 

9.73 ±0.15 

9.75 

± 

0.32 

SDSS J095744.46-I-330820.7 

4.23 

99.00 

47.17 ±0.02 

46.67 

±0.01 

48.05 ± 

-1.00 

9.19 

± 

0.67 

99.00 

± 

99.00 

10.44 ±0.23 

99.00 

± 

99.00 

QUEST J101046.3-013104.1 

5.09 

99.00 

99.00 ±99.00 

47.70 

±0.01 

48.02 ± 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

SDSS J101336.37-I-561536.4 

3.63 

3.64 

46.95 ±0.03 

46.99 

±0.09 

45.58 ± 

0.14 

6.90 

± 

0.91 

6.65 

± 

1.88 

10.06 ±0.24 

10.32 

!± 

:0.31 

SDSS J102622.89-I-471907.0 

4.94 

99.00 

99.00 ±99.00 

46.64 

±0.01 

47.97 ± 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

SDSS J103221.11-1-092749.0 

4.00 

4.02 

47.02 ±0.01 

46.55 

±0.02 

45.31 ± 

0.33 

7.58 

± 

0.43 

99.00 

± 

99.00 

10.18±0.13 

99.00 

± 

99.00 

SDSS J103242.71-I-433605.4 

3.46 

99.00 

46.50 ±0.07 

46.40 

±0.01 

47.62 ± 

-1.00 

5.97 

± 

0.19 

99.00 

± 

99.00 

9.68 ±0.15 

99.00 

± 

99.00 

SDSS J104351.204-650647.6 

4.47 

99.00 

46.85 ±0.05 

46.46 

±0.01 

47.91 ± 

-1.00 

4.07 

± 

1.14 

99.00 

± 

99.00 

9.52 ±0.32 

99.00 

± 

99.00 

SDSS J104433.04-012502.2 

5.80 

99.00 

99.00 ±99.00 

46.51 

±0.03 

48.35 ± 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

SDSS J104437.06-I-650645.0 

3.65 

99.00 

46.88 ±0.03 

46.60 

±0.01 

47.81 ± 

-1.00 

4.25 

± 

0.23 

99.00 

± 

99.00 

9.57 ±0.20 

99.00 

± 

99.00 

SDSS J104845.05-I-463718.3 

6.19 

99.00 

99.00 ±99.00 

46.55 

±0.03 

48.50 ± 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

SDSS J105123.03-I-354534.3 

4.91 

4.91 

46.86 ±0.05 

46.70 

±0.01 

45.52 ± 

0.18 

3.35 

± 

0.17 

3.94 

±- 

-1.00 

9.35 ±0.20 

9.75 

±- 

-1.00 

SDSS J105756.28-I-455553.0 

4.14 

4.13 

47.43 ±0.02 

47.24 

±0.10 

45.90 ± 

0.09 

6.93 

± 

0.25 

4.22 

± 

0.70 

10.33 ±0.22 

lO.lC 

l± 

:0.23 

SDSS J110657.83-I-081643.3 

4.27 

99.00 

46.79 ±0.02 

46.45 

±0.03 

47.95 ±- 

-1.00 

4.85 

± 

0.35 

99.00 

± 

99.00 

9.65 ±0.15 

99.00 

± 

99.00 

SDSS J111348.13-I-540940.0 

3.78 

3.77 

99.00 ±99.00 

46.57 

±0.01 

45.12 ± 

0.30 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

BRI 1114-0822 

4.49 

99.00 

46.58 ±0.03 

46.19 

±0.02 

48.07 ±- 

-1.00 

6.06 

± 

0.12 

99.00 

± 

99.00 

9.74 ±0.20 

99.00 

± 

99.00 

SDSS J113002.35-I-115438.3 

3.40 

99.00 

46.91 ±0.02 

46.49 

±0.02 

47.64 ± 

-1.00 

6.25 

± 

0.08 

99.00 

± 

99.00 

9.94 ±0.14 

99.00 

± 

99.00 

SDSS J113142.31-I-384854.6 

4.13 

99.00 

99.00 ±99.00 

46.74 

±0.01 

47.94 ± 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

SDSS J113246.50-I-120901.6 

5.17 

99.00 

99.00 ±99.00 

46.52 

±0.03 

48.10 ± 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

SDSS J113307.63-I-522835.5 

3.74 

3.74 

46.71 ±0.04 

46.64 

±0.01 

45.33 ± 

0.20 

7.14 

± 

0.37 

2.50 

±- 

-1.00 

9.95 ±0.21 

9.30 

±- 

-1.00 

SDSS J113418.96-I-574204.6 

3.52 

99.00 

46.98 ±0.03 

46.61 

±0.01 

47.67 ± 

-1.00 

7.90 

± 

0.84 

99.00 

± 

99.00 

10.20 ±0.23 

99.00 

± 

99.00 

SDSS J114514.18-I-394715.9 

4.06 

4.03 

47.02 ±0.02 

46.81 

±0.01 

45.27 ± 

0.32 

7.51 

± 

0.34 

99.00 

± 

99.00 

10.17±0.15 

99.00 

± 

99.00 

SDSS J114657.79-I-403708.6 

5.01 

99.00 

46.54 ±0.08 

46.57 

±0.01 

48.07 ±- 

-1.00 

4.43 

± 

0.71 

99.00 

± 

99.00 

9.43 ±0.25 

99.00 

± 

99.00 

SDSS J114816.64-I-525150.2 

6.42 

99.00 

99.00 ±99.00 

46.70 

±0.02 

48.19 ± 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

SDSS J115757.96-I-485655.7 

4.25 

99.00 

46.72 ±0.05 

46.47 

±0.01 

47.93 ± 

-1.00 

4.72 

± 

0.30 

99.00 

± 

99.00 

9.58 ±0.21 

99.00 

± 

99.00 

SDSS J115935.63-I-042420.0 

3.45 

99.00 

46.53 ±0.04 

46.12 

±0.03 

47.62 ±. 

-1.00 

4.10 

± 

0.36 

99.00 

± 

99.00 

9.35 ±0.16 

99.00 

± 

99.00 
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TABLE 7 — Continued 


Name 

(1) 

-^ref 

(2) 

^Ha 

(3) 

log Li350 

(4) 

log ^'5100 

(5) 

log Lho: 

(6) 

FWHM3,c,v 

(7) 

FWHM3,hc< 

(8) 

log Mbh.Civ 

(9) 

log MbH.Hq: 
(10) 

SDSS J120110.31+211758.5 

4.58 

99.00 

47.00 

±0.04 

46.63 

±0.01 

48.03 ± 

-1.00 

12.01 

l± 

:0.48 

99.00 

± 

99.00 

10.60 ±0.22 

99.00 

± 

99.00 

SDSS J120207.78+323538.8 

5.29 

5.27 

47.05 

±0.03 

46.54 

±0.02 

45.35 ± 

0.35 

lo.oe 

l± 

:1.83 

99.00 

± 

99.00 

10.46 ±0.28 

99.00 

± 

99.00 

SDSS J120441.71-002149.5 

5.03 

99.00 

46.48 

±0.07 

46.45 

±0.03 

47.99 ± 

-1.00 

5.71 

± 

1.75 

99.00 

± 

99.00 

9.63 ±0.34 

99.00 

± 

99.00 

SDSS J120447.15-I-330938.7 

3.62 

3.64 

46.93 

±0.15 

46.97 

±0.01 

45.65 ± 

0.13 

8.08 

± 

0.28 

7.82 

± 

2.01 

10.18±0.17 

10.52 

!± 

:0.28 

BR 1202-0725 

4.69 

4.69 

47.32 

±0.02 

46.96 

±0.01 

45.56 ± 

0.07 

9.44 

± 

0.45 

4.78 

± 

0.65 

10.55 ±0.23 

lo.oe 

l± 

:0.20 

SDSS J120934.54-I-553745.4 

3.57 

3.58 

47.13 

±0.02 

46.96 

±0.07 

45.38 ± 

0.14 

3.77 

± 

0.99 

2.50 

±- 

-1.00 

9.60 ±0.31 

9.47 

±- 

-1.00 

SDSS J122738.30-I-572748.9 

3.99 

99.00 

99.00: 

±99.00 

46.79 

±0.01 

47.93 ± 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

SDSS J123239.29-I-525250.9 

4.29 

4.27 

47.08 

±0.02 

46.74 

±0.01 

45.38 ± 

0.28 

7.71 

± 

0.19 

99.00 

± 

99.00 

10.23 ±0.21 

99.00 

± 

99.00 

SDSS J124306.55-I-530522.0 

3.57 

3.55 

46.97 

±0.03 

46.53 

±0.01 

44.95 ± 

0.35 

6.65 

± 

1.02 

99.00 

± 

99.00 

10.04 ±0.25 

99.00 

± 

99.00 

SDSS J125051.93-1-313021.9 

6.13 

99.00 

99.00: 

±99.00 

46.46 

±0.01 

48.45 ±- 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

SDSS J125051.93-1-313021.9 

6.13 

6.09 

99.00: 

±99.00 

46.46 

±0.01 

44.84 ± 

0.34 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

SDSS J130002.16-I-011823.0 

4.61 

99.00 

46.88 

±0.04 

46.67 

±0.02 

48.02 ±- 

-1.00 

5.44 

± 

0.87 

99.00 

± 

99.00 

9.80 ±0.21 

99.00 

± 

99.00 

SDSS J130348.94-I-002010.4 

3.65 

3.62 

46.74 

±0.01 

46.73 

±0.01 

45.57 ± 

0.11 

2.67 

± 

0.12 

6.98 

± 

1.56 

9.54 ±0.08 

10.29 

l± 

:0.26 

SDSS J130608.26-I-035626.3 

5.99 

99.00 

99.00: 

±99.00 

46.36 

±0.11 

48.47 ±- 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

SDSS J131914.20-I-520200.0 

3.90 

3.93 

47.22 

±0.02 

46.85 

±0.01 

45.41 ± 

0.25 

4.65 

± 

0.20 

99.00 

± 

99.00 

9.84 ±0.21 

99.00 

± 

99.00 

SDSS J132420.83-I-422554.6 

4.04 

4.01 

46.85 

±0.04 

46.65 

±0.01 

45.47 ± 

0.16 

6.76 

± 

1.00 

3.33 

± 

1.22 

9.99 ±0.25 

9.57 

± 

0.37 

SDSS J132423.26-I-623342.1 

3.63 

99.00 

46.48 

±0.06 

46.35 

±0.01 

47.68 ± 

-1.00 

5.69 

± 

0.13 

99.00 

± 

99.00 

9.63 ±0.20 

99.00 

± 

99.00 

SDSS J132603.00-I-295758.1 

3.77 

99.00 

46.93 

±0.03 

46.23 

±0.02 

47.72 ±- 

-1.00 

7.30 

± 

0.69 

99.00 

± 

99.00 

10.10 ±0.23 

99.00 

± 

99.00 

SDSS J133223.26-I-503431.3 

3.81 

99.00 

47.01 

±0.02 

46.56 

±0.01 

47.72 ±- 

-1.00 

6.89 

± 

1.06 

99.00 

± 

99.00 

10.09 ±0.25 

99.00 

± 

99.00 

SDSS J133412.56-I-122020.7 

5.13 

99.00 

46.68 

±0.06 

46.36 

±0.03 

48.16 ± 

-1.00 

4.38 

± 

0.97 

99.00 

± 

99.00 

9.49 ±0.28 

99.00 

± 

99.00 

SDSS J133448.70-I-521317.9 

3.61 

3.60 

46.71 

±0.05 

46.51 

±0.01 

45.17± 

0.25 

3.68 

± 

0.25 

99.00 

± 

99.00 

9.35 ±0.20 

99.00 

± 

99.00 

SDSS J133529.45-I-410125.9 

4.26 

4.29 

47.21 

±0.01 

46.82 

±0.01 

45.78 ± 

0.15 

6.45 

± 

1.72 

3.89 

± 

1.21 

10.14 ±0.29 

9.80 

± 

0.32 

BRI 1335-0417 

4.40 

99.00 

46.81 

±0.02 

46.44 

±0.02 

48.08 ±- 

-1.00 

10.33 

;± 

:0.29 

99.00 

± 

99.00 

10.35 ±0.22 

99.00 

± 

99.00 

SDSS J134015.03-I-392630.7 

5.03 

5.05 

46.68 

±0.06 

46.48 

±0.02 

45.15 ± 

0.30 

8.94 

± 

0.53 

99.00 

± 

99.00 

10.15 ±0.22 

99.00 

± 

99.00 

SDSS J134040.24-I-281328.1 

5.34 

5.38 

99.00: 

±99.00 

46.54 

±0.02 

45.56 ± 

0.23 

99.00 

± 

99.00 

2.94 

± 

1.45 

99.00 ±99.00 

9.39 

± 

0.48 

SDSS J134743.29-I-495621.3 

4.51 

4.53 

47.39 

±0.02 

46.97 

±0.10 

45.51 ± 

0.19 

7.13 

± 

1.63 

6.78 

± 

2.34 

10.33 ±0.30 

10.39 

l± 

:0.36 

BRI 1346-0322 

3.99 

3.99 

46.64 

±0.03 

46.82 

±0.01 

45.26 ± 

0.16 

4.77 

± 

0.08 

99.00 

± 

99.00 

9.55 ±0.19 

99.00 

± 

99.00 

SDSS J140132.76-I-411150.4 

4.02 

99.00 

46.48 

±0.04 

46.51 

±0.01 

47.88 ±- 

-1.00 

5.49 

± 

1.62 

99.00 

± 

99.00 

9.59 ±0.30 

99.00 

± 

99.00 

SDSS J140146.53-I-024434.7 

4.44 

99.00 

47.01 

±0.02 

46.65 

±0.02 

48.00 ± 

-1.00 

6.37 

± 

2.69 

99.00 

± 

99.00 

10.00 ±0.33 

99.00 

± 

99.00 

SDSS J140248.08-I-014634.0 

4.16 

99.00 

47.09 

±0.03 

46.62 

±0.01 

48.07 ±- 

-1.00 

8.11 

± 

1.47 

99.00 

± 

99.00 

10.28 ±0.27 

99.00 

± 

99.00 

SDSS J141831.70-1-444937.5 

4.31 

4.33 

47.04 

±0.02 

46.54 

±0.01 

45.32 ± 

0.29 

6.10 

± 

0.39 

99.00 

± 

99.00 

10.00 ±0.16 

99.00 

± 

99.00 

SDSS J142144.98-I-351315.4 

4.56 

99.00 

46.87 

±0.05 

46.54 

±0.01 

47.88 ±- 

-1.00 

6.98 

± 

1.33 

99.00 

± 

99.00 

10.03 ±0.27 

99.00 

± 

99.00 

SDSS J142242.13-I-461310.2 

3.73 

99.00 

46.98 

±0.03 

46.50 

±0.00 

47.79 ± 

-1.00 

7.76 

± 

0.57 

99.00 

± 

99.00 

10.18 ±0.22 

99.00 

± 

99.00 

SDSS J142243.02-I-441721.4 

3.55 

3.61 

47.05 

±0.01 

47.18 

±0.05 

45.71 ± 

0.07 

11.87 ± 

:0.53 

6.11 

± 

0.84 

10.61 ±0.16 

10.41 

.± 

:0.21 

SDSS J143835.95-I-431459.2 

4.61 

4.67 

47.35 

±0.01 

47.14 

±0.05 

45.39 ± 

0.20 

9.84 

± 

0.42 

4.92 

± 

1.89 

10.60 ±0.14 

10.19 

l± 

:0.39 

SDSS J144144.76-I-472003.2 

3.63 

3.62 

46.71 

±0.06 

46.56 

±0.01 

45.24 ± 

0.18 

3.31 

± 

0.23 

2.84 

± 

1.32 

9.26 ±0.20 

9.37 

± 

0.45 

SDSS J144213.09-I-391856.0 

3.63 

99.00 

47.01 

±0.02 

46.43 

±0.01 

47.58 ±- 

-1.00 

7.53 

± 

0.76 

99.00 

± 

99.00 

10.17±0.17 

99.00 

± 

99.00 

SDSS J144340.70-I-585653.2 

4.28 

99.00 

47.24 

±0.02 

46.75 

±0.10 

47.96 ± 

-1.00 

9.19 

± 

0.23 

99.00 

± 

99.00 

10.48 ±0.22 

99.00 

± 

99.00 

SDSS J144350.66-I-362315.1 

5.27 

99.00 

46.78 

±0.03 

46.58 

±0.01 

48.08 ±- 

-1.00 

9.46 

± 

1.49 

99.00 

± 

99.00 

10.25 ±0.21 

99.00 

± 

99.00 

SDSS J144542.75-I-490248.9 

3.88 

3.87 

47.30 

±0.02 

47.12 

±0.01 

45.94 ± 

0.07 

3.14 

± 

0.12 

6.57 

± 

0.97 

9.52 ±0.20 

10.44 

1 ± 

:0.21 

SDSS J144733.46-I-465723.6 

3.59 

99.00 

46.59 

±0.06 

46.29 

±0.01 

47.58 ±- 

-1.00 

4.79 

± 

0.90 

99.00 

± 

99.00 

9.53 ±0.26 

99.00 

± 

99.00 

SDSS J145408.95-I-511443.7 

3.64 

3.61 

47.22 

±0.02 

47.08 

±0.08 

45.46 ± 

0.16 

4.64 

± 

0.27 

4.68 

± 

1.52 

9.84 ±0.21 

10.11 

.± 

:0.34 

SDSS J150620.48-I-460642.4 

3.50 

3.52 

46.84 

±0.02 

46.38 

±0.01 

45.21 ± 

0.22 

7.21 

± 

0.21 

6.98 

± 

3.02 

10.03 ±0.15 

10.11 

.± 

:0.42 

SDSS J150654.54-I-522004.6 

4.07 

99.00 

47.11 

±0.02 

46.92 

±0.01 

47.16 ± 

-1.00 

9.24 

± 

1.43 

99.00 

± 

99.00 

10.41 ±0.21 

99.00 

± 

99.00 

SDSS J151035.29-I-514841.0 

5.03 

5.02 

99.00: 

±99.00 

46.47 

±0.01 

45.27 ± 

0.35 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

SDSS J151442.74-I-532412.0 

3.57 

99.00 

46.65 

±0.06 

46.42 

±0.01 

47.58 ±- 

-1.00 

5.77 

± 

0.11 

99.00 

± 

99.00 

9.73 ±0.20 

99.00 

± 

99.00 

SDSS J152034.53-I-383906.5 

3.40 

99.00 

46.47 

±0.05 

46.11 

±0.01 

47.46 ±- 

-1.00 

6.92 

± 

0.46 

99.00 

± 

99.00 

9.79 ±0.15 

99.00 

± 

99.00 

SDSS J152413.35-I-430537.4 

3.93 

3.95 

46.83 

±0.02 

46.58 

±0.01 

45.23 ± 

0.29 

5.07 

± 

0.89 

99.00 

± 

99.00 

9.71 ±0.21 

99.00 

± 

99.00 

SDSS J153650.25-I-500810.3 

4.93 

4.92 

47.01 

±0.07 

46.84 

±0.01 

45.49 ± 

0.17 

2.65 

± 

0.65 

3.49 

± 

1.25 

9.21 ±0.30 

9.71 

± 

0.37 

SDSS J153725.35-014650.3 

3.45 

99.00 

46.55 

±0.06 

46.04 

±0.04 

47.56 ± 

-1.00 

5.62 

± 

0.15 

99.00 

± 

99.00 

9.65 ±0.20 

99.00 

± 

99.00 

SDSS J153825.74-I-420933.3 

3.83 

99.00 

46.82 

±0.02 

46.52 

±0.01 

47.79 ±- 

-1.00 

7.78 

± 

0.60 

99.00 

± 

99.00 

10.10±0.16 

99.00 

± 

99.00 

SDSS J154340.38-I-341744.4 

4.41 

99.00 

47.15 

±0.02 

47.01 

±0.01 

48.02 ± 

-1.00 

5.16 

± 

0.15 

99.00 

± 

99.00 

9.90 ±0.15 

99.00 

± 

99.00 

SDSS J154914.45-I-364822.0 

3.58 

3.60 

46.59 

±0.03 

46.26 

±0.01 

44.88 ± 

0.35 

6.96 

± 

0.51 

99.00 

± 

99.00 

9.87±0.16 

99.00 

± 

99.00 

SDSS J160254.18-I-422822.9 

6.07 

99.00 

99.00: 

±99.00 

46.65 

±0.01 

48.47 ±- 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

SDSS J161140.13-I-273029.6 

3.33 

99.00 

46.70 

±0.04 

46.30 

±0.02 

47.58 ±- 

-1.00 

4.52 

± 

0.38 

99.00 

± 

99.00 

9.53 ±0.16 

99.00 

± 

99.00 

SDSS J161425.13-I-464028.9 

5.31 

99.00 

99.00 

±99.00 

46.59 

±0.01 

48.07 ±- 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

SDSS J162100.70-I-515544.8 

5.59 

5.59 

99.00 

±99.00 

46.72 

±0.01 

45.20 ± 

0.18 

99.00 

± 

99.00 

3.09 

± 

1.17 

99.00 ±99.00 

9.53 

± 

0.38 

SDSS J162331.81-1-311200.5 

6.22 

6.23 

99.00 

±99.00 

46.39 

±0.04 

44.83 ± 

0.21 

99.00 

± 

99.00 

4.37 

±- 

-1.00 

99.00 ±99.00 

9.68 

±- 

-1.00 

SDSS J162520.31-I-225832.9 

3.77 

3.78 

47.11 

±0.01 

46.66 

±0.06 

45.44 ± 

0.17 

5.02 

± 

0.18 

5.37 

± 

1.80 

9.85 ±0.15 

10.01 

.± 

:0.34 

SDSS J162623.38-I-484136.4 

4.89 

99.00 

47.09 

±0.05 

46.69 

±0.01 

47.92 ± 

-1.00 

2.70 

± 

0.20 

99.00 

± 

99.00 

9.27±0.21 

99.00 

± 

99.00 

SDSS J162943.43-I-391211.4 

3.91 

99.00 

46.71 

±0.03 

46.32 

±0.01 

47.93 ± 

-1.00 

7.19 

± 

0.44 

99.00 

± 

99.00 

9.97±0.16 

99.00 

± 

99.00 

SDSS J163636.92-I-315717.0 

4.56 

4.57 

46.95 

±0.03 

46.55 

±0.03 

45.57 ± 

0.23 

4.90 

± 

0.85 

6.66 

± 

2.87 

9.74 ±0.21 

10.15 

.± 

:0.42 

SDSS J163847.42-I-232716.4 

3.82 

3.84 

46.88 

±0.06 

46.95 

±0.01 

45.28 ± 

0.27 

5.09 

± 

0.45 

99.00 

± 

99.00 

9.74 ±0.17 

99.00 

± 

99.00 

SDSS J163950.52-I-434003.7 

3.99 

99.00 

47.28 

±0.01 

46.95 

±0.01 

47.87 ±- 

-1.00 

11.53 

;± 

:0.74 

99.00 

± 

99.00 

10.71 ±0.17 

99.00 

± 

99.00 

SDSS J164248.71-I-240303.3 

3.48 

3.48 

46.96 

±0.02 

46.41 

±0.01 

45.11 ± 

0.22 

6.53 

± 

0.45 

2.50 

±- 

-1.00 

10.01 ±0.16 

9.17 

±- 

-1.00 

SDSS J165354.61-I-405402.1 

4.98 

4.97 

47.05 

±0.01 

46.71 

±0.01 

45.45 ± 

0.08 

6.60 

± 

0.34 

4.36 

± 

0.66 

10.07 ±0.13 

9.85 

± 

0.20 

SDSS J165436.85-I-222733.7 

4.70 

4.68 

47.09 

±0.02 

46.81 

±0.01 

45.48 ± 

0.24 

6.51 

± 

0.61 

99.00 

± 

99.00 

10.08 ±0.17 

99.00 

± 

99.00 

SDSS J165902.11-I-270935.1 

5.31 

99.00 

99.00 

±99.00 

46.58 

±0.02 

48.04 ± 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

SDSS J172100.76-I-601721.0 

5.80 

99.00 

99.00 

±99.00 

46.75 

±0.01 

47.81 ± 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

SDSS J173744.87-I-582829.5 

4.92 

99.00 

99.00 

±99.00 

46.42 

±0.01 

48.07 ± 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

PSS J1745+6848 

4.13 

99.00 

99.00 

±99.00 

46.31 

±0.01 

47.99 ± 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

RX J1759.4+6638 

4.32 

99.00 

99.00 

±99.00 

46.22 

±0.01 

47.96 ± 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

PSS J1802+5616 

4.16 

99.00 

99.00 

±99.00 

46.10 

±0.01 

48.17 ± 

-1.00 

99.00 

± 

99.00 

99.00 

± 

99.00 

99.00 ±99.00 

99.00 

± 

99.00 

BR 2212-1626 

3.99 

3.98 

46.92 

±0.03 

47.03 

±0.01 

45.93 ± 

0.10 

2.31 

± 

0.10 

6.73 

± 

1.26 

9.04 ±0.19 

10.42 

!± 

:0.23 

BRI 2235-0301 

4.25 

99.00 

47.02 

±0.02 

47.12 

±0.01 

47.92 ± 

-1.00 

10.21 

i± 

:0.39 

99.00 

± 

99.00 

10.46 ±0.22 

99.00 

± 

99.00 

SDSS J223841.81-000105.2 

3.49 

3.49 

46.86 

±0.02 

46.55 

±0.02 

44.97 ± 

0.33 

6.10 

± 

0.36 

99.00 

± 

99.00 

9.90 ±0.15 

99.00 

± 

99.00 
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TABLE 7 — Continued 


Name 

-^ref 

^Ha 

log Li 350 

log ^5100 

log Lho: 

FWHM3,c,v 

FWHM 3 ,hc< 

log Mbh,Civ 

log Mbh,H a 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 


(7) 

(8) 

(9) 

(10) 

BR 2248-1242 

4.16 

4.14 

46.66 ±0.07 

46.81 ±0.01 

45.82 ± 

0.10 

1.58 ±0.07 

4.74 ±0.96 

8.55 ±0.19 

9.98 ±0.24 

BR J2328-4513 

4.36 

99.00 

99.00 ±99.00 

46.41 ±0.01 

47.92 ±- 

-1.00 

99.00 ±99.00 

99.00 ±99.00 

99.00 ±99.00 

99.00 ±99.00 

SDSS J233330.17-1-152538.7 

3.68 

3.68 

46.61 ±0.06 

46.31 ±0.02 

45.19 ± 

0.26 

6.07 ±0.88 

99.00 ±99.00 

9.75 ±0.24 

99.00 ±99.00 

BR J2349-3712 

4.21 

99.00 

46.86 ±0.02 

46.53 ±0.01 

47.95 ±- 

-1.00 

5.73 ±0.37 

99.00 ±99.00 

9.84 ±0.21 

99.00 ±99.00 


Note. — Catalog of the properties derived for the AKARI quasars, sorted by right ascension. Column 1: Target name; Column 2: Redshift 
from references; Column 3: Redshift measured from Ho:; Column 4: 1350 A luminosity and its uncertainty; Column 5: 5100 A luminosity and its 
uncertainty; Column 6: Ha luminosity and its uncertainty; Column 7: FWHM of the C IV line and its uncertainty; Column 8: FWHM of the Ha 
line and its uncertainty; Column 9: Mbh from the C IV line and its uncertainty; Column 10: Mbh from the Ha line and its uncertainty. The units 
for L, FWHM, and Mbh are ergss”^, 1000kms~^, and Mq. Columns 9 and 10 are from Equations (10) and (7), respectively. Empty parameters 
are entered as 99 and upper limits are given with errors of -1. 



